For  Reference 


NOT  TO  BE  TAKEN  FROM  THIS  ROOM 


(3x  mm 

WliffiiE 


THE  UNIVERSITY  OF  ALBERTA 


Release  Fora 


NAME  OF  AUTHOR  Praveer  Asthana _ 

TITLE  OF  THESIS  An  Inverse  Method  in  One  Dimension 

DEGREE  FOR  WHICH  THESIS  WAS  PRESENTED  M.Sc. _ 

YEAR  THIS  DEGREE  GRANTED  1982 


Permission  is  hereby  granted  to  THE  UNIVERSITY  OF 

ALBERTA  LIBRARY  to  reproduce  single  copies  of  this  thesis 

•  * 

and  to  lend  or  sell  such  copies  for  private  scholarly  or 
scientific  research  purposes  only. 

The  author  reserves  other  publication  rights,  and 
neither  the  thesis  nor  extensive  extracts  from  it  may 
be  printed  or  otherwise  reproduced  without  the  author’s 
written  permission. 


THE  UNIVERSITY  OF  ALBERTA 


AN  INVERSE  METHOD  IN  ONE  DIMENSION 


BY 


PRAVEER  ASTHANA 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 
OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  PHYSICS 


EDMONTON,  ALBERTA 


FALL,  1982 


THE  UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES  AND  RESEARCH 


The  undersigned  certify  that  they  have  read,  and  recommend  to 
the  Faculty  of  Graduate  Studies  and  Research  for  acceptance,  a  thesis 
entitled  AN  INVERSE  METHOD  IN  ONE  DIMENSION  submitted  by  PRAVEER  ASTHANA 
in  partial  fulfilment  of  the  requirements  for  the  degree  of  Master  of 


Science. 


ABSTRACT 


Given  a  sequence  of  bound-state  energies,  the  inverse  Gel'fand- 
Levitan  theory  for  the  one-dimensional  Schrodinger  equation  enables  us 
to  uniquely  construct  a  symmetric  reflectionless  potential  which  sup¬ 
ports  those  bound-states.  A  question  of  interest  is  that  if  the  bound- 
state  energies  belong  to  a  symmetric,  but  not  necessarily  reflectionless, 
potential  then  how  closely  does  the  reflectionless  potential  simulate 
the  actual  one? 

This  method,  if  reliable,  gives  us  a  means  of  approximately 
constructing  a  symmetric  'non-confining  well'  in  one  dimension  from  its 
bound-state  spectrum.  With  the  inclusion  of  an  additional  parameter, 

'E  ' ,  this  technique  can  also  be  used  for  'locally'  approximating  a 
symmetric  'confining'  potential  in  one  dimension. 

In  order  to  test  the  reliability  of  this  technique,  we  apply  it 
to  some  known  symmetric  potentials  (confining  and  non-confining)  in 
one  dimension.  We  also  calculate  the  reflection  coefficients  of  the 
actual  potentials  (of  the  'local'  non-confining  counterparts  in  the 
case  of  the  confining  ones)  as  a  function  of  momentum  and  look  for  some 
'favourable'  (or  'unfavourable')  systematics  towards  this  end. 

In  conclusion,  the  reflectionless  potential  simulates  the  rough 
qualitative  features  of  the  actual  potential  quite  well.  Within  this 
perspective,  the  technique  survives  the  tests. 
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1.  INTRODUCTION 


The  'November  revolution',  even  in  the  winter  of  1974,  had 
made  things  hum  in  High-Energy  Physics.  A  new  particle  ip/J  (with  a 
mass  of  3.096  GeV)  was  discovered  at  the  Brookhaven  National  Laboratory 
in  the  study  of  the  reaction  [1] 

p  +  Be  (e+e  )  +  anything 

and  almost  simultaneously  at  SLAC  [2]  in  the  reaction 
e+e  -+  hadrons  . 

Its  companion  particle  ip'  (3.686  GeV)  was  almost  immediately  spotted 
at  SLAC  [3].  Amidst  growing  theoretical  expectations,  it  confirmed 
the  need  for  the  introduction  of  a  fourth  'c(aharn) -quark'  in  the  des¬ 
cription  of  hadronic  matter.  It  was  accepted,  almost  immediately,  that 
ip  and  ipr  were  different  energy  states  of  a  system  made  up  of  a  c-quark 
and  a  c-antiquark.  Since  their  discovery  in  November  1974,  a  rather 
satisfactory  spectrum  of  cc  bound-states  has  been  unfolded  [4-7] .  A 
recently  reported  spectrum  is  being  tabulated  in  Table  1.1  [4,5]. 

The  discovery  of  cc  states  stimualted  the  quest  for  such  bound- 
systems  of  quarks  with  still  newer  'flavours'.  In  the  spring  of  1977, 
barely  two  and  a  half  years  after  the  entry  of  ip  to  the  hadronic  world, 
another  particle  'upsilon  (T) '  was  discovered  at  Fermilab  [8,9]  in  the 
reaction 

p  +  N  -*  (y+  +  y  )  +  anything. 

Almost  immediately,  it  was  found  to  be  accompanied  by  another  partner 
T'  [10] .  It  has  been  a  field  of  intense  experimentation  since  then. 
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Table  1.1.  Charmonium  States 


State 

Mass  (MeV) 

n 

c 

2978 

±  4 

3096.9  ±  0.1 

3412 

±  5 

*1 

3508.4  ±  4 

x2 

3553.9  ±  4 

K 

3592 

±  5 

r 

3686 

±  0.2 

Table  1.2.  Bottomonium  States 


State 

Mass  (MeV) 

T 

9459  ±  10 

T'  -  T 

560  ±  3 

rplt  _  rp 

889  ±  4 

rplll  _  rp 

1114  ±  5 

Note:  The  numbers  tabulated  above  are  being  quoted  from  a  very  recent 
paper.  The  reference  is:  K.J.  Miller  and  M.G.  Olsson 3  Phys. 
Rev.  D25  (1982)  2383. 


So  far  four  companions  of  T  have  been  detected  experimentally  [5,11] 
as  shown  in  Table  1.2.  In  much  the  same  way  as  ip  and  ip',  these  par¬ 
ticles  have  been  interpreted  as  different  energy  states  of  a  bound- 
system  of  a  new  quark,  ’b (bottom) '  and  its  antiquark,  b. 

As  the  cc  and  bb  systems  resemble  positronium,  in  the  sense 
of  being  a  bound  state  of  a  particle  and  its  antiparticle,  they  are 
known  as  the  "charmonium1  and  ’bottomonium ?  respectively.  The  generic 
name  ’ quarkonium'  has  been  assigned  to  a  system  of  this  kind. 

It  is  enough  of  a  testimony  of  the  significance  of  the  cc-  and 
bb-systems  that  they  have  attracted  special  attention  in  a  medley  of 
hundreds  of  mesons.  What  is  so  fascinating  about  them?  Let  us  con¬ 
sider  the  cc-system  for  instance.  It  is  a  very  heavy  meson  built  up 
of  two  very  heavy  quarks.  The  binding  energy  available  to  the  quarks 
is  very  small  compared  to  their  mass.  So,  it  is  a  system  which  can  be 
analyzed  using  the  tools  of  the  usual  non-relativistic  quantum  mecha¬ 
nics,  unlike  the  case  of  the  systems  made  up  of  lighter  quarks.  The 
situation  is  still  better  and  more  promising  in  bottomonium  as  it  is 
still  heavier  in  mass.  To  highlight  this  latitude  in  dealing  with 
these  particles,  people  have  called  these  systems  as  the  "hydrogen- 
atoms  of  strong-interaction  physics  [12]".  It  is  worthwhile  to 
acknowledge  another  remark  [7]:  "The  discovery  of  charm  marks  the 
beginning  of  quark  chemistry". 

At  the  core  of  this  enthusiasm  is  the  hope  to  unravel  the  nature 
of  interquark-potential  using  simple  non-relativistic  quantum  techni¬ 
ques  on  these  systems.  Since  the  discovery  of  charmonium  in  1974,  a 


lot  of  effort  has  gone  into  this  direction.  These  efforts  were  further 
invigorated  by  the  discovery  of  bottomonium  in  1977.  It  still  remains 
the  motivation  for  a  very  active  phenomenology. 

The  majority  of  the  work  done  in  this  direction  can  be  broadly 
classified  into  two  categories.  In  the  first  case,  guided  by  theoreti¬ 
cal  suggestions,  one  assumes  a  specific  form  of  spherically  symmetric 
static  potential.  This  consists  of  some  free  parameters.  One  then 
tries  to  calculate  the  level  schemes  and  other  known  relevant  quanti¬ 
ties  using  the  Schrodinger  equation.  Fitting  these  to  the  known  spec¬ 
trum  and  quantitites,  one  fixes  the  parameters.  There  is  a  whole  variety 
of  such  potentials  which  do  the  job  [6,13-21].  With  certain  forms  of 
potential,  e.g.  the  power  law  potentials,  the  Schrodinger  equation  dic¬ 
tates  simple  scaling  laws  for  many  quantities  of  interest  like  level 
spacings,  decay  widths,  etc.  This  is  the  topic  of  interest  in  the 
second  category.  By  comparing  the  predicted  and  observed  systematics 
in  these  quantities  one  can  fix  the  details  of  the  potential  or,  at 
least,  establish  some  bounds  on  its  characteristics.  In  addition  to 
these,  many  general  quantum  mechanical  and  semi-classical  results  have 
also  been  used  to  extract  information  regarding  the  inter quark-potential 
from  the  available  experimental  data  [6,20-35]. 

What  one  has  to  realize  is  that  in  all  such  calculations,  one 
has  to  assume  some  kind  of  a  potential  -  no  matter  how  general  it  may 
be  -  to  start  with.  It  will  be,  thus,  highly  desirable  and  satisfying 
to  have  a  method  which  can  give  us  the  ’confining’  interquark  potential 
straight  from  the  measured  quantities.  The  measured  quantities  are 
the  bound-state  energies  and  the  square  of  the  bound-state  wavefunctions 
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at  origin  which  are  known  from  the  leptonic  decay  widths  [6,20],  This 
defines  an  inverse  qucan.tiMi-meohanical  problem. 

Such  aspirations  are  not  very  recent  in  the  history  of  physics. 
Determination  of  a  molecular  potential  from  its  bound-state  spectrum 
using  WKB  inversion  techniques  dates  back  to  the  late  1920’s  and  early 
’30’s  [36-44].  An  application  of  these  techniques  to  the  present  prob¬ 
lem  will  require  that  the  bound-state  energy  E  be  considered  as  a 

'  n 

continuous  function  of  the  principal  quantum  number  ’n1  and  the  deri¬ 
vative  (dE^/dn)  be  known  as  an  input  [20,33,45].  This  information  is 
not  available  to  us.  Moreover,  the  known  energy  levels  are  so  few  in 
number  that  a  polynomial  fit  to  E^  as  a  function  of  n  by  the  standard 
numerical  methods  will  not  be  highly  dependable  [33]  .  In  the  present 
context,  therefore,  one  needs  a  different  method  to  achieve  this  goal. 


Quigg,  Rosner  and  Thacker  (QRT)  came  forward  with  an  ingenious 

3 

suggestion  [45-47] .  They  considered  the  n  levels  of  the  quarkonium. 
The  well-known  correspondence  between  the  radial  Schrodinger  equation 
with  Z = 0  and  the  one-dimensional  Schrodinger  equation,  then,  tells  us 
that  these  are  nothing  but  the  odd-parity  levels  of  a  symmetric  poten¬ 


tial  in  one  dimension.  The  square  of  the  n  wavefunction  at  the 

i  1 2 

origin,  1^(0)  |  ,  can  also  be  related  to  the  odd-parity  wavefunction  in 
one  dimension. ^2n> the  well-known  result  [20,46,47] 


'U°>l2  -  ^  ir2n(°>l2 


n 


(1.1) 


where  the  prime  denotes  the  derivative  with  respect  to  ’x’.  QRT’s  idea 
was  to  look  for  an  inverse  method  for  the  one-dimensional  Schrodinger 
equation  which  can  uniquely  reconstruct  a  symmetric  potential  from  its 


odd-parity  bound-state  energies  and  the  slopes  of  the  odd-parity  wave 
functions  at  the  origin. 


In  fact,  they  dug  out  of  the  literature,  a  modest  and  practic¬ 
able  way.  It  has  been  known  for  a  long  time  that  a  reflectionless 
potential  can  be  uniquely  constructed  in  terms  of  2N  parameters,  where 

N  is  the  number  of  bound  states  it  can  support  [45,47,49-56],  N  of 

2 

these  parameters  are  the  bound-state  energies  =  -K^.  If  one  demands, 
in  addition,  that  the  potential  be  symmetric,  i.e. 

V (x)  =  V(-x)  (1.2) 

then  the  other  N  parameters  also  get  related  to  the  N  bound-state 
energies  in  the  following  way  [45,47] 


n 

m^n 


K  +  K 

m  n 


K  -  K 

m  n 


(1.3) 


The  symmetric  reflectionless  potential  is  given  by  [45,47,49] 


V(x)  =  -2  ^ -j  {£n  det  (I+A)} 
dx 


(1.4) 


I  is  the  N*N  unit  matrix  and  A  is  the  matrix  given  by 
c  ,1/2  ,1/2  e(Kn+Kv)x 

A  =  An  K  VtTT  ' 

n  v 


(1.5) 


The  normalized  bound-state  wavefunctions  of  this  potential  are  given  by 


H,  (X)  =  _ i _  det(l_+A)_ 

U  a1^2  ^KnX  det  (I+A) 


(n) 


(1.6) 


n 


where  (l+A)^n^  is  the  matrix  (I+A)  with  its  n^  column  differentiated 


once  with  respect  to  'x'.  Thus,  we  can  reconstruct  a  symmetric  reflec- 
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tionless  potential  in  one  dimension  only  from  its  bound-state  spectrum. 


However,  this  method  in  its  present  form,  is  still  not  suf¬ 
ficient  to  serve  the  purpose  in  the  context  of  the  quarkonium.  In  a 
quarkonium  problem,  one  has  to  approximate  a  'confining’  potential 
for  which  the  whole  edifice  of  scattering  theory  breaks  down.  It  has 
an  infinite  number  of  bound  states  and  no  scattering  states  at  all. 

The  last  of  QRT's  suggestions  was  to  remove  this  deadlock.  According 
to  them,  one  should  try  to  construct  a  N  bound-state  approximation 
V^(x)  to  the  confining  potential  V(x),  where  [V^Cx)  -E  ]  is  a  reflec¬ 
tionless  potential  supporting  N  bound-states  at  the  energies 
2  th 

-K  =  -E  +E  .  Here  E  is  the  n  bound  state  energy  of  the  confining 
n  o  n  n 

potential  and  Eq  is  a  parameter  which  defines  the  zero  of  energy. 

There  is  no  physical  reasoning  which  will  enable  us  to  fix  this  para¬ 
meter  uniquely.  It,  thus,  gives  rise  to  an  2? -  ambiguity .  However, 
it  can  lie  in  the  range 


<  E 


N+l 


(1.7) 


Increasing  N  in  steps  from  N  =  l,  one  can  then  successively  build  up 
'local  approximations'  to  the  confining  potential.  QRT  carried  out 
such  a  local  construction  of  the  harmonic  oscillator,  the  linear 
potential  and  the  infinite  square  well.  Their  numerical  experiments 
by  varying  Eq  in  its  allowed  interval  showed  that  a  consistent  choice 
of 


E 

o 


+Vi> 


(1.8) 


gives  extremely  good  agreement  between  the  confining  potential  and  the 
local  reflectionless  approximation  V^(x) .  In  fact,  the  approximation 


seemed  to  converge  to  the  actual  potential  ’locally'  as  N  became 
larger  [45,47] . 

In  the  application  of  this  technique  to  the  actual  three- 
dimensional  inverse  problem,  the  odd  parity  levels  of  the 

corresponding  symmetric  potential  in  one  dimension  are  given  by  the 
energy  levels  of  the  three-dimensional  potential.  A  knowledge  of  the 
slopes  of  the  odd-parity  wave  functions  at  the  origin  then  yields 
simultaneous  equations  for  the  'unphysical'  even-parity  levels 
via  eqn.  (1.6).  The  solutions  can  be  obtained  by  numerical  methods 
[46,47,48],  Thus  one  has  all  the  bound-states  of  a  symmetric  poten¬ 
tial  in  one-dimension.  Using  the  inverse  technique  outlined  above, 
one  can  then  approximate  it  by  a  reflectionless  potential.  Such  an 
application,  in  the  context  of  quarkonium,  was  undertaken  in  a  prac¬ 
tical  way  by  QRT .  They  were  able  to  obtain  the  charmonium  and  the 
bottomonium  potentials  [46]  .  They  also  used  such  a  technique  to  in¬ 
vestigate  the  ' flavour- independence '  of  the  interquark  potential  [48]. 

Quarkonia  are  not  the  only  systems  where  such  an  inversion 
might  be  desired.  The  most  attractive  feature  of  this  technique  is 
that  it  gives  us  a  direct,  though  approximate,  information  regarding 
the  potential  from  a  few  bound-state  energies.  Such  a  technique  can, 
thus,  have  a  much  wider  applicability.  These  handsome  prospects  were 
precisely  the  motivation  for  further  studying  this  technique  in  this 
thesis.  The  task  undertaken  in  this  thesis  is  to  test  the  reliability 
of  this  technique  by  applying  it  to  a  larger  variety  of  known  symmetric 
potentials  in  one  dimension.  In  Chapter  2  we  develop  the  formal  as¬ 
pects  of  this  technique.  Chapter  3  is  then  devoted  to  the  applications. 


' 


It  was  realized  that  confining  potentials  are  not  very  clean 
cases  to  serve  as  tests  of  this  technique.  There  is  an  inherent,  but 
unavoidable,  artificiality  involved  in  translating  the  pure  bound- 
state  problem  of  a  confining  potential  to  a  scattering  problem.  This 
technique,  in  its  own  right,  is  applicable  only  to  ’local*  scattering 
(non-confining)  potentials.  Moreover,  not  all  the  potentials  of 
interest,  in  general,  are  confining.  So,  one  of  the  tasks  undertaken 
in  this  thesis  is  to  estimate  some  local  symmetric  potentials  in  one 
dimension  by  the  corresponding  reflectionless  potentials  using  this 
technique.  Since  we  are,  in  general,  approximating  a  reflecting  po¬ 
tential  by  a  reflectionless  one,  it  would  be  interesting  to  calculate 
its  reflection  coefficient  and  look  for  some  systematics  in  its 
behaviour  as  a  function  of  momentum  which  suggest  a  close  agreement 
between  the  actual  and  the  estimated  potentials.  Section  3.1  is 
devoted  to  this  task. 

In  Section  3.2  we  estimate  some  confining  potentials  using  the 
QRT  prescription. 

Section  3 .3  is  devoted  to  a  discussion  concerning  the  conver¬ 
gence  of  the  reflectionless  approximation  to  the  actual  reflecting 
potential.  Our  aim  is  to  try  to  unravel  various  factors  which  affect 
the  convergence,  both  ’favourably'  and  ’unfavourably’. 

Chapter  4  concludes  this  work  with  a  critical  discussion  of 
the  results  obtained;  it  weighs  the  enthusiasm  regarding  the  prospects 
of  this  technique  against  the  skepticism. 
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2.  THE  TECHNIQUE  AND  ITS  NUMERICAL  TEST 


In  this  chapter,  we  first  briefly  review  the  Gel ’f and-Levitan 
formalism  of  inverse  scattering  in  one  dimension.  Then  we  derive  the 
technique  for  approximating  a  symmetric  potential  by  a  symmetric  re¬ 
flectionless  potential  having  the  same  bound-state  spectrum.  Finally, 
we  test  this  technique  by  applying  it  to  a  known  symmetric  reflection¬ 
less  potential  in  one  dimension. 

2.1  (DIRECT)  SCATTERING  IN  ONE  DIMENSION 

2.1.1  One-dimensional  Schrodinger  Equation  and  the  Fundamental  Solutions. 

We  start  with  the  one -dimensional  Schrodinger  equation  Cft = 2m = 1) 
2 

^  +  [k^  -  V(x)]u(x,k)  =  0  .  (2.1) 

dx 

2 

Here  u(x,k)  is  the  wavef unction,  V(x)  is  the  potential  and  k  (=E)  is 
the  energy  of  the  particle.  The  variable  x  is  defined  on  the  infinite 
interval,  -°°  <  x  <  oo. 


If  the  potential  is  localized,  i.e.,  if 

V(x)  -►  0  as  [ x I  oo  (2.2) 


then  the  solutions  u(x,k)  of  the  Schrodinger  equation  will  reduce  to  a 

+  ilex 

linear  combination  of  plane  waves  e  at  x  = ±°o.  It  is,  therefore, 
convenient  to  express  all  solutions  of  the  Schrodinger  equation  as  linear 
combinations  of  the  two  solutions  f^(x,k)  and  f?(x,k)  which  exhibit  the 
following  large  x  behaviour 


f1(x,k)  =  e 


ikx 


f2(x,k)  =  e 


-ikx 


as  x  +°o' 


as  x  -*■  — °° . 


(2.3) 
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These  solutions  are  commonly  referred  to  as  the  fundamental  solutions 
of  the  Schrodinger  equation. 

2.1.2  Integral  Equations  for  the  Fundamental  Solutions. 

It  is  obvious  by  looking  at  eqn.  (2.1)  that  explicit  forms  of 
f^(x,k)  and  f^Cxjk)  cannot  be  obtained  unless  we  are  given  the  potential 
V(x).  However,  if  we  treat  k  as  a  complex  variable,  their  analytic 
properties  in  the  complex  k-plane  can  be  inferred  from  their  integral 
f  orms . 


Using  the  Green  function  for  eqn.  (2.1),  we  can  immediately  get 
the  following  integral  equations  for  f^(x,k)  and  f^(x,k) 


f1(x,k)  =  e 


ikx  1 


roo 


dx'  sin  k(x-x')  V(x’)  f^(x',k)  (2.4) 


f2(x,k)  =e-ikx+i 


X 


dx’  sin  k(x-x')  V(x?)  f^(x'  ,k)  .  (2.5) 


If  the  potential  V(x)  is  real,  then  we  get  the  following  two  properties 


of  f^(x,k)  and  f^(x,k)  from  their  integral  equations 

JU  JL 

f1(x,k)  =  f ^(x ,-k  ) 


f2(x,k)  =  f2(x,-k  ) 


(2.6) 


If,  in  addition,  k  is  also  real,  then  we  get  from  eqn.  (2.6) 


f^(x,k)  =  f ^(x, -k) 


* 


f  2  (x  ,k)  =  f  2  (x  ,-k) , 


(2.7) 


Let  us  now  try  to  solve  these  integral  equations  by  iteration. 
In  eqn.  (2.4),  for  example,  we  first  put  e  for  f^(x,k)  and  so  on. 


. 
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The  resulting  integrals  would  then  be  found  to  converge  for  lmk>0.  It  is 
also  well-known  that  the  resulting  series  for  Volterra  integral  equations 
such  as  eqn.  (2.4)andeqn.  (2.5)  arealways  uniformly  convergent  with  respect 
to  k  [57.58],  Thus,  we  conclude  that  f^(x,k)  is  an  analytic  function  in  the 
upper  half  of  the  complex  k-plane.  Inexactly  the  same  way,  f2(x,k)  will 
also  be  an  analytic  function  in  the  upper  half  of  the  complex  k-plane. 

2.1.3  Definitions  of  the  Reflection  and  the  Transmission  Coefficients. 

We,  now,  define  the  Wronskian  of  two  solutions  of  the  Schrodinger 
equation  as 

W[u1;u2]  S  u^u2  -  u.^  .  (2.8) 

This  definition  differs  from  the  conventional  one  by  a  minus  sign.  Here 
the  prime  means  a  dif f erentation  with  respect  to  x. 

The  Schrodinger  equation  (2.1)  does  not  contain  the  first  deri¬ 
vative  of  u(x,k).  So,  the  Wronskian  of  any  two  linearly  independent 
solutions  will  be  independent  of  x  but  it  may  depend  on  k  [59] .  We  can 
exploit  this  property  in  calculating  the  following  useful  Wronskian 
relations 

W[f 1(x,k) ; f 1(x,-k) ]  =  W[f1(x,k);f1(x,-k)]x=oo  =  2ik  (2.9) 

W[f 2(x,k) ; f 2(x,-k) ]  =  W[f2(x,k) ;f2(x,-k)]x=_oo  =  -2ik.  (2.10) 

These  Wronskian  relations  depend,  obviously,  only  on  k.  Hence, 

{f ^(x ,k) ,f ^(x,-k) }  and  (f2(x,k) ,f2(x,-k) }  form  two  sets  of  linearly  in¬ 
dependent  solutions  to  the  Schrodinger  equation  on  the  infinite  interval. 


We  can,  thus,  write 


■ 
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f 2 (x ,k)  =  c11(k) 

f1(x,k)  +  c12(k)  f1(x,-k) 

(2.11) 

f1(x,k)  =  c21(k) 

f2(x,-k)  +  c22 (k)  f 2 (x ,k) 

(2.12) 

These  two  equations,  as  we  discuss  below,  describe  two  complementary 
scattering  situations  in  one  dimension. 


Using  the  asymptotic  forms,  eqns.  (2.3),  of  f^(x,k)  and  f^(x,k) 

we  conclude  from  eqn.  (2.11)  that  it  represents  a  wave  e  at  x  =  -°° 

_  ikx 

and  a  linear  combination  c^(k)e  +  C-^^e  at  x  =  +°°.  Therefore, 
this  represents  a  scattering  situation  in  which  a  wave  of  amplitude 
C12^10  is  incident  from  the  right  on  the  potential  V(x).  It  is  re¬ 
flected  with  an  amplitude  C-Q^k)  and  transmitted  with  an  amplitude  unity 

___  % 

to  x  =  -°°.  This  has  been  drawn  schematically  in  Figure  2.1. 


Similarly,  eqn.  (2.12)  represents  a  scattering  situation  in  which 
a  wave  of  amplitude  c2]/k)  is  incident  from  the  left.  A  part  c22^)  is 
reflected  and  a  wave  of  unit  amplitude  is  transmitted  to  x  =  +°°.  Figure 
2.2  represents  this  scattering  from  the  left,  schematically. 


The  familiar  reflection  coefficient ,  R(k) ,  and  the  transmission 
coefficient ,  T(k) ,  for  waves  of  unit  incident  amplitude  can  now  be  de¬ 
fined  as  follows 


cn(k) 

c12tk) 


Vk)  ■ 


(2.13) 


c.„00 


‘  Vk)  ■  z^oo 


(2.14) 


where  the  subscripts  R  and  L  refer  to  the  scattering  of  waves  incident 
from  the  right  and  the  left  respectively. 


Fig.  2.1.  A  wave  incident  from  the  right. 


c21(Weikx 

7 

Scattering 

Region 

ikx 

e 

 c22(k)e-ikx 

/ 

Fig.  2.2.  A  wave  incident  from  the  left. 


2.1.4  Relationships  between  these  Coefficients. 


A  number  of  relationships  exist  between  the  various  coefficients 
defined  in  the  previous  section.  It  is  particularly  illuminating  to 
establish  them  for  some  of  those  relations  are  statements  of  important 
physical  realities. 

Let  us  substitute  f^(x,k)  from  eqn.  (2.12)  into  eqn.  (2.11)  and 
equate  the  coefficients  of  f^(x,k)  and  f2(x,-k)  on  both  the  sides.  We 
finally  get  the  following  two  relations 


ll(k) 

c22°° 

+  c12(k) 

c2i(-k) 

u(k) 

c2i(k) 

+  C12^k^ 

c22(-k)  . 

(2.15) 

Similarly,  if  we  substitute  f^(x,k)  from  eqn.  (2.11)  into  eqn.  (2.12) 
and  equate  the  coefficients  of  f^(x,k)  and  f^(x,-k)  on  both  the  sides, 
we  get 

1  =  c21(k)  c12(-k)  +  c22^  cn^ 

0  =  c21(k)  e^C-lc)  +  c22^  C12^k^  *  (2-16) 

Equations  (2.9)  to  (2.12)  enable  us  to  establish  the  following  expres¬ 
sions  for  c . . ' s 
iJ 

cll(k)  =  2^-  W[f2(x,k) ;f1(x,-k) ]  (2.17) 

c22(k)  =  2lk  w[f2^x,-k);f1(x,k)]  (2.18) 

c12 (k)  =  C21^k^  =  Ilk  wtfi^x»k^ ;f2^x»k^ •  (2.19) 


Combining  the  expressions  for  T^(k)  and  T^(k)  with  eqn.  (2.19),  we 


.  Ik  .  >4' 


immediately  conclude  that 


TL(k)  =  TR(k)  *  (2.20) 

This  is  commonly  referred  to  as  the  reprooity  relation. 

In  the  event  when  V(x)  and  k  are  real  i.e.  when  eqn.  (2.7)  holds 
good,  one  can  show  that 

cl2(-k)  =  ci2 (k)  (2.21) 

and 

cn(k)  =  _c22  (k)  =  ~c22  (“k)  .  (2.22) 

Using  eqn.  (2.21)  and  eqn.  (2.22)  in  eqn.  (2.15)  and  eqn.  (2.16),  we  get 

lcl2(k)|2  =  1  +  |cu(k)|2  =  1  +  |c22(k)|2  .  (2.23) 

By  the  definitions  of  the  reflection  and  transmission  coefficients,  eqn. 
(2.23)  implies 

1  =  |  T(k)  | 2  +  |  Rr  (k)  | 2  =  j  T(k)  |  2  +  iR^k)!2  (2.24) 

where,  in  view  of  the  reciprocity  relation  (2.20),  we  have  used  T(k) 

for  both  the  left  and  right  transmission  coefficients.  Eqn.  (2.24)  is 
the  statement  of  the  law  of  conservation  of  energy . 

If  one  starts  with  the  definitions  (2.13)  and  (2.14)  and  uses 
the  relationships  (2.21)  and  (2.22),  one  gets  the  following  equations 

R^k)  T(-k)  +  R^-k)  T(k)  =  0  .  (2.25) 

Eqn.  (2.25)  is  commonly  referred  to  as  the  phase  law. 


It  is  also  a  simple  matter  to  verify,  using  eqn.  (2.21)  and 


- 


1  ■ 
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eqn.  (2.22)  that 

R^(-k)  =  R^k) 

R^(-k)  =  P^Ck)  (2.26) 


2.1.5.  A  Bridge  between  these  Coefficients  and  the  S-Matrix  in  One 
Dimension. 

We  have  already  remarked  that  remote  from  a  localized  potential 

any  solution,  u(x,k) ,  of  the  Schrodinger  equation  in  one  dimension,  can 

+  ikx 

be  written  as  a  linear  combination  of  the  plane  waves  e  .  Therefore, 
we  have  [60] 

lim  (A^e^^  +  B+e  -  u(x,k)}  =  0  (2.27) 

x-*+°°  "r 

and 

lim  {A  e^*  +  B  e  -  u(x,k)}  =  0  (2.28) 

x+-°° 


We  define  an  incoming  wave  as  a  right  moving  wave  at  x  = -00  and  a  left 

moving  wave  at  x  =  4^°.  Similarly  an  outgoing  wave  is  a  right  moving 

wave  at  x  = +°°  or  left  moving  wave  at  x  = -00.  So,  the  incoming  parts  of 

•j  —  iloc 

u(x,k)  are  A  e  and  B+e  whereas  the  outgoing  parts  of  u(x,k)  are 

.  ikx  ,  n  -ikx 
A+e  and  B  e 


We  can  further  define  a  two-dimensional  vector  space  of  column 
vectors.  The  first  component  of  each  vector  represents  the  amplitude 
of  the  right  moving  part  and  the  second  represents  the  left  moving  part . 
So ,  the  incoming  wave  vector  is 

/  N 

A 


- 


whereas  the  outgoing  wave  vector  is 


r  \ 

A+ 


B 

V.  _ > 


Right  from  the  definitions  of  the  reflection  and  the  transmission  co¬ 
efficients,  it  is  obvious  that  a  relationship  must  exist  between  these 
two  vectors.  To  establish  that,  note  that  the  incoming  wave  vector  has 
a  wave  of  amplitude  A  incident  from  x  = -°°.  This  will  get  reflected 
with  an  amplitude  R^(k)A  and  transmitted  with  an  amplitude  T^(k)A_. 

So,  this  will  give  rise  to  an  outgoing  wave  R^(k)A  at  x  = -00  and  another 
T  (k)A  at  x=-H».  Similarly,  we  have  a  wave  of  amplitude  B  incident 

J_j  “  i 

from  the  right.  This  will  give  rise  to  an  outgoing  wave  of  amplitude 
R^(k)B+  to  the  right  and  another  one  of  amplitude  TR(k)B+  to  the  left. 
So,  the  outgoing  wave  vector  has  a  component  [T  (k)A  +  R  (k)B  ]  moving 
towards  x  =  +°°  and  another  component  [R_^(k)A_  +  T^(k)B+]  moving  towards 
x  =  — 00 .  In  other  words. 


«J+ 

TL(k)A_  +  RR(k)B+' 

B 

,RL(k)A_  + 

i.e. 


f  \ 

A+ 


=  S 


r  \ 

A 


(2.29) 


where  S  is  a  2x2  matrix  given  by 


S  = 


'V» 

Vk) 


RR(k) 

Vk) 


(2.30) 


We  can  alternatively  write  eqn.  (2.29)  as 
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u^(x,k)  =  S  ^(x,k)  (2.31) 

where  u^+^(x,k)  and  \x,k)  stand  for  the  outgoing  and  incoming  parts 
of  the  wave  function  respectively.  The  matrix  S  is  the  so-called  S- 

matrix . 

This  establishes  the  correspondence  between  the  reflection  and 
the  transmission  coefficients  and  the  more  common  description  of  the 
scattering  theory  in  terms  of  the  S-matrix.  It  is  obvious  that  the 
relations  (2.24)  and  (2.25),  derived  in  the  previous  section,  confirm 
the  desired  unitarity  of  the  S-matrix. 

2.1.6  Bound-State  Poles  of  the  Transmission  Coefficient. 

The  potential  V(x)  can,  in  general,  support  bound  states.  Where 
does  this  characteristic  of  the  potential  get  reflected?  This  is, 
clearly,  a  vital  information  for  the  problem  we  are  concerned  with. 

And,  in  this  section,  we  explore  this  question. 

Let  us  start  by  investigating  the  analytic  structure  of  the 
transmission  coefficient  in  the  complex  k-plane.  We  first  prove  that 
the  poles  of  the  transmission  coefficient  in  the  upper  half-plane  lie 
on  the  positive  imaginary  axis.  By  eqn.  (2.13)  and  eqn.  (2.14),  it  is 
obvious  that  we  are  interested  in  the  zeros  of  c^C^)*  By  eqn.  (2.24), 
which  holds  only  for  real  k,  it  is  evident  that  c12^  can  never  vanish 
on  the  real  k-axis,  i.e. 

c  2(k)  ^  0  if  Imk  =  0  . 

So,  the  zero  of  must  lie  off  the  real  axis. 


(2.32) 


V  J 


. 


Let  be  the  location  of  one  of  the  poles.  The  Schrodinger 
equation  gives 

uM  +  k2u  =  V(x)u  (2.33) 

o 

and  its  complex  conjugate  dictates 

u  "  +  k  u  =  V(x)u  (2.34) 

o 

where  we  have  assumed  that  V(x)  is  real.  u(x,k)  may  be  either  of  the 

* 

two  fundamental  solutions  f^(x,k)  or  f2(x,k)  and  u  (x,k)  is  the  cor- 

■k 

responding  complex  conjugate.  Multiplying  eqn.  (2.33)  by  u  and  eqn. 
(2.34)  by  u,  and  subtracting 

(u  u"  -  u  u)  +  (k  -  k  *)u  u  =  0  .  (2.35) 

o  o 


Integrating  over  all  x 


f°° 

dx 

*  — oo 


2  ,  *  * 
u |  =  I (u  u '  -  uu  ' ) 


00 

—CO 


(2.36) 


The  right  hand  side  of  eqn.  (2.36)  is  either  of  the  two  Wronskians  in 
eqn.  (2.10),  both  of  which  are  independent  of  x.  So,  it  is  zero. 
Therefore , 


(k2  -  k*2> 
o  o 


dx 


(2.37) 


If  we  write 


k  =  k  +  ik  . 
o  or  oi 


then  eqn.  (2.37)  becomes 


k  k  . 
or  oi 


dx 

— oo 

1  u  [ 

0 


Thus,  if  k  .  i  0,  then  k 
’  oi  or 


(2.38) 


(2.39) 


=0,  because  the  integral  is  a  positive  and 
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finite  quantity  (we  prove  this  claim  below)  when  lmk>0. 


If  c^Ck  )  t*ien  from  ecln.  (2.11)  and  eqn.  (2.12) 

f  (x,k  )  =  c  (k  )  f  (x,k  ) 

2  o  11  o  1  o 

f  (x,k  )  =  c  (k  )  f  (x,k  ) 

I  o  22  o  2  o 

Hence,  c  (k  )  and  cn_(k  )  are  related  by 
11  o  22  o 

C-i  -i  (k  )  =  - - - r-  . 

II  o  c00 (k  ) 

22  o 


(2.40) 


(2.41) 


So,  the  two  solutions  are  linearly  dependent.  From  the  asymptotic 

forms  of  f,(x,k)  and  f0(x,k)  we  infer  that  if  Imk  >0,  then  f,(x,k  ) 

1  2  o  1  o 

exponentially  decreases  at  x  =  +00  while  the  solution  f_(x,k  )  behaves 

2  o 

likewise  at  x  = -00.  Thus,  for  k  =  i«  ,  we  have  a  solution  which  is 

o  o 

quadrat ically  integrable  on  the  entire  axis  [55]  .  The  characteristic 
exponential  decay  shows  that  f^(x,iKQ),  or  equivalently  f2(x,iKQ),  cor¬ 
responds  to  a  bound-state  solution. 


2.1.7.  The  Bound-State  Wavefunctions  and  their  Normalizations. 

With  the  results  of  the  preceding  section  at  our  disposal,  we 

know  that  if  the  potential  V(x)  supports  N  bound  states  at  k =  iK^ 

(n = 1,2 , . . . ,N) ,  then  the  transmission  coefficient  has  N  poles  on  the 

positive  imaginary  axis  at  k^  =  iK^  (n  =  1,2, . . .  ,N).  The  nt^1  bound-state 

wavefunction  is  given  by  f,(x,iK  )  or  f^(x,iK  )  with  the  relation 

°  y  1  n  2  n 

f2(x,iKn)  =  c11(iKn)  f^XjiiO  (2.42) 

where  c..,(iK  )  and  c__(iK  )  are  related  by 
11  n  22  n 

c_.(iK  )  =  - j-. — s’ 

11  n  c00(ik  ) 

22  n 


(2.43) 


- 
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Let  us  now  try  to  derive  expressions  for  the  normalization  constants 

of  f,(x,iK  )  and  f0(x,iK  )  . 

-L  n  2  n 


As  the  first  step,  we  calculate  c10(k  )  defined  by 

12  n  J 


S9(k  >  E 

12  n 


dci2(k) 

dk 


(2.44) 


k=kn 


Using  eqn.  (2.19)  9  we  get 


^12 (k)  = 


- ~~2  W[f1(x,k)  ;f2(x,k)  ]  + (w[  f  f  2  ]  +W[f1;f2]} 


2ik 


By  virtue  of  eqn.  (2.42),  the  first  term  on  the  right-hand  side  vanishes 


at  k  =  k  and  the  second  reduces  to 
n 


C12^kn^  2ik  ^Cll^kn^  W^l;fl^  +C22^kn^  ^2,45') 
n 

We  have  to  now  evaluate  the  Wronskians  on  the  right-hand  side  of  eqn. 
(2.45).  Towards  this  end,  we  start  by  writing  the  Schrodinger  equa¬ 
tions  corresponding  to  f^(x,k)  and  f^(x,kn) 


fl^(x,k)  +  k  f^(x,k)  =  V (x)  f1(x,k) 


(2.46) 


f  "(x ,k  )  +  k  f  (x,k  )  =  V (x)  f .  (x ,k  ) 
1  *  n  n  1  n  I  n 


(2.47) 


Multiplying  eqn.  (2.46)  by  f  (x,k  )  and  eqn.  (2.47)  by  f^Xjk)  and  then 
subtracting  one  equation  from  the  other 

^  [f^(x,kn)f1(x,k)  -£|(x,k)f1(s,kn)]  -  (k2 -k^)f1Cx,k)f1(x,kn)=0 

(2.48) 


Differentiating  eqn.  (2.48)  with  respect  to  k  and  then  putting  k  =  k^, 

-j-  W[f ,(x,k  );f  (x,k  )]  =  2k  [f  (x,k  )]2  . 
dx  l  n  l  n  nl  n 


(2.49) 


. 
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Integrating  eqn.  (2.49)  from  x  to  +°°,  we  get 


W[f  (x',k  );f  (x\k  )] 

-i-  n  i  n 


00  ,00 

=  2k 


X 


X 


dx' [f  (x' ,k  ) ]  .  (2 .50) 

i  n 


It  is  obvious  from  eqn.  (2.3)  that  for  lmkn>0,  f^(x,kn)-*-0  as  x^+°°, 
Hence,  the  Wronskian  in  eqn.  (2.50)  will  vanish  at  x  = +°°  yielding 


W[f, (x  k  );f  (x,k  )]  =  -2k 
1  n  1  n  n 


x 


dx? [f1(xl ,kn) ] 


(2.51) 


Following  exactly  the  same  reasoning  for  f2(x,k^),  we  get 


W[f  (x,k  ) ; f  (x,k  )]  =  2k 
z  n  z  n  n 


rx 


_oo 


dx* [f  (x* ,k  ) ] 
2  n 


(2.52) 


Substituting  eqn.  (2.51)  and  eqn.  (2.52)  into  eqn.  (2.45),  we  get 


roo 


c  (k  )  =  — 
12  n  i 


_oo 


dx  f . (x,k  )f  (x,k  ) 
1  n  2  n 


(2.53) 


As  an  immediate  corollary  to  eqn.  (2.53) ,  we  can  see  that  the 

zeros  of  (k)  for  a  real  potential  can  only  be  simple.  First  of  all, 

from  the  integral  equations  for  f^(x,k)  and  f2(x,k)  we  see  that  for 

real  potential  and  Imk  >0,  f. (x,k  )  and  f„(x,k  )  are  real.  Now,  eqn. 
r  n  1  n  2  n 

(2.53)  can  be  written  as. 


c  (k  )  =  -ic  (k  ) 
12  n  11  n 


fOO 

dx[f1(x,kn) ]2  . 

— OO 


Since  fn (x,k  )  is  real,  the  integrand  on  the  right  hand  side  is  positive 

and  the  integral  is  non-zero.  As  we  will  see  below,  the  normalization 

constant  of  f, (x,k  )  is  proportional  to  ^  (k  ).  So,  cni(k  )  ^0.  Hence, 
1  *  n  11  n  11  n 

c, _ (k  )  ^0.  This  implies  that  the  zeros  of  c._(k),  or  equivalently,  the 
12  n  lz 

poles  of  the  transmission  coefficient,  have  to  be  simple  for  a  real  po¬ 


tential. 


4 
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Eqn.  (2.53)  leads  us  to  the  residue  of  the  transmission  coefficient  at 

k  =  k  ,  viz . 
n* 


iy  =  Residue  of 
n 


cl200 


c.,(k  ) 
12  n 


— oo 


dx  f  (x,k  )f  (x,k  ) 
1  n  2  n 

(2.54) 


Using  eqn.  (2.42)  and  eqn.  (2.43),  we  obtain  from  eqn.  (2.54),  the  fol¬ 
lowing  two  relations 


f  oo 


Y. 


n 


/•OO 


— OO 


dx  f  (x ,k  )f  (x,k  )  =  1  = 
1  n  2  *  n 


—OO 


dx  y  c  (k  ) [f  (x,k  ) ] 
n  11  n  In 


(2.55) 


and 


Y. 


OO 


rco 


—oo 


dx  f-(x,k  )f  (x,k  )  =  1  = 
1  n  2  n 


dx  y  c  9(k  ) [f  (x,k  )] 
n  22  n  2  n 


(2.56) 


1/2 

It  is  evident  from  eqn.  (2.55)  and  eqn.  (2.56)  that  [Ync-Q(kn^  and 

1/2 

[y  c__(k  )]  are  the  normalization  constants  for  the  bound-state  wave- 
n  zz  n 

functions  f^(x,k^)  and  f2(x,kn),  respectively. 

Also,  since  f1_(x,kn)  and  f2^x,kn^  are  rea1,  ir  follows  from 
eqn.  (2.54)  that  Y0  is  real.  We  can  write  the  normalization  constants 

Xs 


as : 


c  (k  ) 

V  E  Vn(V  =  -1  dTrnrr 

1Z  n 


LJ 


_oo 


dx[f1(x,kn) ] 


-1 


(2.58) 


and 


C22(kn) 

Yn  E  V22(kn)  '  ^  d^TFT 


r°° 


lj  _oo 


dx[f  (x,k  )] 
Z  n 


-1 


.  (2.59) 


As  f _ (x ,k  )  and  f.(x,k  )  are  real,  A  and  A  are  also  real  and  positive. 
1  *  n  2  n  Kn  uu 

Using  the  results  obtained  above,  one  can  immediately  prove  the  following 


[55] 


• 

t  S' 


* 


. 
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(2.60) 


Proof  of  eqn.  (2.60):  We  have  remarked  above  that  c_.(k  )  ^  0  due  to 

11  n 

the  normalization  condition.  So: 


L.H.S.  =  Y 


2 

n 


=  y  c  (k  )  * 
n  11  n 


Y. 


n 


cTTTk) 

11  n 


=  [y  c  (k  )][Y  c  0(k  )]  ,  from  eqn.  (2.40) 
n  11  n  n  ZZ  n 


R.H.S. 


QED. 


2.2  INVERSE  SCATTERING  IN  ONE  DIMENSION 

2.2.1.  The  Prelude:  A  Representation  of  the  Fundamental  Solutions  and 
the  Potential  in  terms  of  a  Function  K(x,y)  Characterizing  the 
Scattering . 

We  have  been  considering  the  time- independent  Schrodinger  equa¬ 
tion  so  far.  The  solutions  of  this  equation  are  monochromatic  waves. 

The  scattering  of  an  incident  pulse  can  then  be  constructed  by  a  Fourier 
synthesis  of  the  results  of  steady-state  calculations.  However,  there 
are  some  advantages  of  considering  the  scattering  of  pulses  directly 
in  the  time  domain  [56,60].  Firstly,  we  are  led  to  a  new  representation 
for  the  fundamental  solutions,  f^(x,k)  and  f2(x,k) ,  in  terms  of  a  func¬ 
tion  K(x,y)  characterizing  the  scattering  and  secondly,  we  are  able  to 
relate  the  scattering  potential  V(x)  to  the  function  K(x,y).  This  is 
the  main  aim  of  this  section. 


The  Fourier  transform  u(x,t)  of  u(x,k)  is  a  transient  wave 


) 


' 


solution  of  the  equation 


9  u(x,t) 
3x2 


1  9  u(x ,t) 

2  2 
c  9t 


-  V(x)  u(x  ,t)  =  0 


(2.61) 


Obviously,  if  the  potential  is  zero  or  if  it  is  localized,  then  remote 
from  the  potential,  eqn.  (2.61),  looks  like 


9  u(x,t)  1  9  u(x,t)  =  n 

2  ~  2  2  -  u 
3x  c  9t 


(2.62) 


An  exact  solution  of  eqn.  (2.62)  is  6(t  ±— )  .  However,  since  there  is 
a  scattering  potential  V(x)  ,  we  would  expect  this  delta  function  pulse  to  be 
partly  scattered  and  thus  leave  some  sort  of  disturbance  behind  a  sharp 
leading  edge.  We,  thus,  write  a  solution  to  eqn.  (2.61)  as 

u2(x,t)  =  6(t+^-)  +  c9(t+^-)  K^^-ct)  .  (2.63) 


Here  0(x)  is  the  unit  step  function 
carries  the  scattering  information, 
has  dimensions) . 


and  K^(x,-ct)  is  the  function  which 
Factor  c  is  just  a  constant  (but 


If  we  take  the  Fourier  transform  of  eqn.  (2.63)  we  get,  with 

ct  =  x  ’  and  k  =  —  , 

c 


U2(x,u)) 


c°° 

,  io)t  .  x 
dt  e  u2(x,t) 

— .OO 


dx»  e_ikx  K^(x,x  ’) ;  x'<x  .  (2.64) 

— OO 

It  can  be  easily  shown  that  U2(x,io)  is  a  solution  to  the  following  equa¬ 
tion 


9  U2  (x  ,co) 
9x2 


+ 


'OJ 


vc 


-  V  (x) 


U2  (x  ,co)  =  0  . 


(2.65) 
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Also,  it  is  obvious  from  eqn.  (2.64)  that 

*i  Icy 

lim  [e  U  (x,co)]  =  1  .  (2.66) 

x+— 00 


Remembering  that  (oj/ c)  =  k,  we  see  that  eqn.  (2.65)  is  nothing  but  the 
time- independent  Schrodinger  equation  (2.1).  We  know  that  f^Xjk)  is  a 
unique  solution  eqn.  (2.65)  or  equivalently  eqn.  (2.1),  with  the  same 
asymptotic  form  as  in  eqn.  (2.66).  Hence,  it  follows  that  I^XjO))  is 
another  representation  of  f^(x,k) 


(x  ,w)  =  f ^  (x  ,k)  =  e  ikx  + 


rx 


—00 


— -j  Ify  ^ 

dx’  K^(x,x?)e  ;  x1  <x  . 

(2.67) 


Similarly,  the  other  solution  to  eqn.  (2.61)  is 

u^(x,t)  =  6(t  -^)  +  c@(t-^-)  KR(x,ct)  .  (2.68) 


Its  Fourier  transform  is  given  by 

ikx  f°°  ikx f 

U1(x,o3)  =  e  +  dx'  K^XjX ' ) e  ; 


x  ’  >  x  .  (2 .69) 


Eqn.  (2.69)  is  again  a  solution  of  the  time-independent  Schrodinger 
equation  with  the  asymptotic  property 

lim  [e~lkx  U1(x,w)]  =  1  .  (2.70) 

X-*co 


Therefore,  U^(x,co)  must  be  another 


ikx 


U1(x,oj)  =  f^(x,k)  =  e  + 


representation  of  f^(x,k) 

r°°  ikxf 

dx’  K  (x,x')e1KX  ;  x'  >x.  (2.71) 
R 
J  x 


These  alternative  representations  of  the  fundamental  solutions  will  be 


extremely  useful  in  our  discussion  of  the  inverse  problem. 


As  we  remarked  above,  all  scattering  information  is  stored  in  the 
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function  K^(x,y)  or  K^(x,y)  for  the  left  and  right  scattering  situations 
respectively.  We  now  show  that  the  scattering  potential  V(x)  is  related 
to  K^(x,x’)  or  K^(x,x’)  by  a  simple  relation.  Substituting  eqn.  (2.63) 
into  eqn.  (2.61),  we  get 

23KL(x,-ct)  ^  3K^(x,-ct) 


-<5(t+-) 

c 


+  c6(t  +— ) 
c 


3x 


3t 


+  V(x) 


32KL(x,-ct)  1  32KL(x,-ct) 


3x 


3t‘ 


-  V(x)KL(x,-ct) 


=  0  . 
(2.72) 


X  X 

Let  us  integrate  eqn.  (2.72)  from  t=- - 8  to  t=-— +e.  As  e+0,  the 

c  c 

second  term  on  the  left  hand  side  makes  no  contribution.  The  delta 
function  in  the  first  term  gives  the  following  non-zero  contribution 

3K^(x,-ct) 


-2 


9K^(x ,-ct) 


3x 


-2 


-c  t=x 


3(-ct) 


+  V(x)  =  0  .  (2.73) 


-ct=x 


And,  obviously  this  relates  V(x)  with  K^(x,x’)  by 
V(x)  =  2  ~  K^(x,x) . 


(2.74) 


Similarly,  if  we  substitute  eqn.  (2.68)  in  eqn.  (2.61)  and  carry 
out  the  same  calculation  as  above ,  we  will  get 


v(x)  =  "2  kr(x’x) 


(2.75) 


Thus,  if  we  have  a  means  of  calculating  K^(x,y)  or  K^(x,y)  from 
the  S-matrix  elements,  we  can  get  the  scattering  potential  via  eqn .  (2.74) 
or  eqn.  (2.75)  respectively.  It  has  been  shown  [55]  that  the  potential 
must  satisfy  the  following  criterion  for  the  validity  of  this  procedure 

[CO 

dx(l  +  |  x  |  )  |  V(x)  |  <  00  (2.76) 

_0O 


‘ 

:  P 
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This  procedure  is  inapplicable  if  the  above  integral  diverges.  So, 
this  excludes  the  cases  where  the  potential  is  either  too  singular  such 
as  V(x)  =6’(x),  or  too  slowly  converging  for  large  values  of  x  [56]. 
And,  if  the  above  criterion  holds  good  then  potentials  obtained  from 
each  of  the  two  equations  (2.74)  and  (2.75)  coincide.  In  fact,  if  we 
are  confident  about  the  validity  of  eqn.  (2.76)  then  we  need  to  know 
only  K^(x,x')  or  K^(x,x')  for  the  inverse  procedure. 

The  final  step  in  establishing  an  inverse  technique  is,  there¬ 
fore,  to  look  for  a  relationship  which  can  give  us  K^(x,y)  or  KR(x,y) 
in  terms  of  the  reflection  and  the  transmission  coefficients.  This  is, 
in  fact,  the  Gel  ' fand-Levitan  Equation. 


2.2.2.  The  Finale;  The  Gel Tf and -Levitan  Equation  for  K(x,y). 

This  integral  equation  can  be  derived  in  a  number  of  ways  [51-53, 
55,56,60].  We  prefer  the  one  using  a  dispersion  technique  [45,50].  We 
start  by  considering  an  analytic  function 


(k) 


f ^(x  ,k)  e 


-ikx 


Imk  >  0 


(2.77a) 


u(x,k)  =  - 


*  *  _  ilcv 

f2(x,k  )e  ;  Imk  <  0  .  (2.77b) 

The  form  in  eqn.  (2.77a)  is  dictated  by  the  fact  that  [l/c^^)  ]f^(x,k) 
gives  the  desired  scattering  boundary  conditions.  As  discussed  in 
Section  (2.1.3),  this  represents  an  incoming  wave  from  the  left  with 
amplitude  unity  and  transmitted  and  reflected  waves  with  amplitudes 
T^(k)  and  R^(k)  respectively. 


Obviously,  by  construction,  u(x,k)  has  poles  on  the  positive 


' 


imaginary  axis  corresponding  to  the  bound— state  poles  of  the  transmission 
coefficient.  Apart  from  these  u(x,k)  has  a  discontinuity  along  the 
real  axis.  The  discontinuity  can  be  easily  calculated  and  it  turns  out 
to  be 


lim[u(x,k+ie)  -u(x,k-ie)]  =  FL  (k)f  (x,k)e 
e-K)  L  1 


—  i  lev 

=  p(x,k).  (2.78) 


The  discontinuity  across  the  cut  along  the  real  axis  is,  therefore, 
proportional  to  the  reflection  coefficient.  In  fact,  this  was  the 
motivation  for  choosing  the  form  in  eqn.  (2.77b). 


From  the  integral  forms  (2.4)  and  (2.5)  for  f^(x,k)  and  f2(x,k), 

we  have 

f^(x,k)  =  eikx  ;  k  =  ±°° 

-ilex 

f 2 (x ,k)  =  e  ;  k  =  ±«  .  (2.79) 


Using  this  in  eqn.  (2.19),  we  get 

cl2(k)  =  1  ;  k  =  ±°° 

Therefore,  it  follows  from  eqn.  (2.77)  that 


u(x,k)  =  1 


;  k  =  ±°° 


(2.80) 


(2.81) 


We  will  thus  construct  u(x,k)  from  its  singularities  using  a  dispersion 
relation  normalized  to  unity  at  a  subtraction  point  at  00 . 


Let  us  consider  u(x,k)  in  the  upper-half  of  the  complex  k-plane, 

We  can  write  the  dispersion  integral 

u(x,k)  _  1  I  u(x,kt)dk' 

[lmk>0]  2t\±  k?-k 


(2.82) 


* 
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where  C  is  the  contour  shown  in  Figure  2.3.  In  the  upper  half-plane 


u(x,k)  has  simple  poles  on  the  positive  imaginary  axis  at  iKn(n=l>2 , . . ,N) 
corresponding  to  the  N  bound  states.  Assuming,  as  justified  above, 
that  u(x,k)  goes  to  zero  at  large  k,  we  can  deform  the  contour  C  in  the 


way  shown  in  Figure  2.4.  We  can,  thus,  write 


u(x,k)  _  1 

[lmk>0]  2iri 


rOO 


— oo 


u(x,k')dk' 

k'-k 


27Ti 


O' 


Y 


u(x  ,k’)dk 1 
k'-k 


(2.83) 


The  contribution  to  the  second  integral  on  the  right-hand  side  comes 

from  the  bound-state  poles.  It  is,  in  fact,  the  sum  of  the  residues  of 

(2.77a)  at  the  bound-state  poles  with  an  overall  minus  sign  coming  from 

the  clockwise  direction  of  the  contour  y.  This  calculation  yields 

K„x 

(2.84) 


i  f  (  i,  in  ju  i  N  f, (x,iK  )e  n 
1  I  u(x,k')dk'  =  y  1  n _ 

27Ti  J  k'-k  L  c_  0  (iK  )  (k-iK  T 

J  y  n=l  12  n  n 


Substituting  eqn.  (2.84)  in  eqn.  (2.83),  we  get 


f  ,  x  N  f  (x ,  iK  )  eKnX  -  r°° 

u(x,k)  =  y  1  n _  1_ 

[ lmk>0]  ^  c10(iK  )  (k-iK  )  2iri 

n=l  12  n  n  J  -°° 


u(x,k')dk' 
k'-k 


(2.85) 


u(x,k)  is,  by  construction,  an  analytic  function  in  lmk<0. 
close  the  contour  C'  as  in  Figure  2.5,  we  have 
1 


0  = 


2TTi 


r 


C' 


u(x,k')dk'  _  1 

k'-k  ”  2ni 


O' 


C  ' 


u(x,k')dk' 

k'-k 


So,  if  we 


(2.86) 


as  the  point  k  is  outside  the  contour.  Since  the  contribution  from  the 
semi-circle  is  zero  eqn.  (2.86)  gives  us 


roo 


o  = 


2iri 


— OO 


u(x,k')dk' 
k'-k 


(2.87) 


In  eqn.  (2.85),  k'  is  at  the  upper  lip  of  the  cut  along  the  real  axis 
whereas  in  eqn.  (2.87)  ,  it  is  at  the  lower  lip  of  the  cut.  If  we 
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Fig.  2.3. 


Fig.  2.4. 


Imk 


Imk 

* 


x 

k 


Fig.  2.5. 
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remember  this  and  subtract  eqn.  (2.87)  from  eqn.  (2.85),  we  have 

K  X 

(  .  N  N  f  (x,iK  ) e  n 

u(x,k)_  y  1  »  n _  _ 

[lmk>0]  L  c10(iK  )  (k-iK  )  2TTi 
n-1  12  n  n 

where  we  have  used  eqn.  (2.78)  for  the  discontinuity  in  u(x,k)  across 
the  cut. 


p(x,k')dk' 

k'-k 


(2.88) 


The  dispersion  relation  in  eqn.  (2.88)  holds  also  for  the  case 
ImkcO.  In  order  to  normalize  u(x,k)  to  unity  at  infinity  we  write 
eqn.  (2.88)  with  one  subtraction  at  infinity 

K  X 


N 


n 


u(x,k)  =  1 +  £  - 


fl(x>iKn)  +  1 


.  c10(iK  )  (k-iK  )  ’’’  2TTi 
n=l  12  n  n 


— CO 


P(Xk'-kdk’-(2-89) 


Using  eqn.  (2.14)  in  eqn.  (2.78),  one  gets 

c22(k,)  -ik’x 

P(x’k,)  =  c^X'k’T  f2(x'k’)e 


Evaluating  the  left  hand  side  of  eqn.  (2.89)  for.lmk<0  using  eqn. 
(2.77b)  and  substituting  the  above  expression  for  p(x,k’)  on  the  right 


hand  side,  we  get 


f2(x,-k)e 


N  f n  (x,iK  )e 

•lkX  =  1  +  l  T  1 


K  X 
n 


L,  c19(i<  Xk-iK  ) 
n=l  12  n  n 


c22Ck'> 


f  00 


27Ti 


dk’ 


’12 


(pyf2(x-k’)e 


-ik'x 


— OO 


k'-k 


(2.90) 


For  Imk  <0,  one  can  write  the  following, 


-ik'x 

e 

k'-k 


-ikx 

-l  e 


x 


— OO 


i(k-k')y 
e  ay 


(2.91) 


and 


' 
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'X  i(k-iK  )y 

e  n  dy  .  (2.92) 

— OO 

•  ~f  Icy 

Multiplying  both  the  sides  of  eqn.  (2.90)  by  e  and  using  eqn.  (2.91) 
and  eqn.  (2.92),  one  gets  after  a  slight  rearrangement, 

f2(x,-k)  =  eikx  + 


A 


KL(x,y)elky  dy 


(2.93) 


k-ix 


=  l  e 


-ikx 


where 


N  f , (x,iK  )  K  y 

Vx’y)  =  iJi^Ten 


2tt 


/•OO 


c„(k  )  , 

dk’  — — -  f  (x  k')e-1^  y 

c10(kf)  V’  ; 

-OO  12 


(2.94) 


Using  the  representation  for  f^(x,k)  derived  in  Section  (2 . 2 . 1) ,  we  have 


-  z  ,  -ik*x 

f  2  (x  ,k’ )  =  e  + 


rx 


— OO 


KL(x,y)e  ik  y  dy  ;  y<x 


and,  hence 


K  X 

f2(x,iiO  =  e  n  + 


x  K  y 

KL(x,y)e  n  dy 


(2.95) 


— OO 


Using  these  two  relations  and  remembering  eqn.  (2.42),  we  get  from  eqn, 
(2.94) 


(x+y )  +  K^x.y)  + 


x 


K  (x,z)£3  (z+y)dz  =  0  ;  y<x  (2.96) 

Li  Li 


where 


^L(z) 


~  dk  C22(k)  -ikz_.  »  c22(iKn)  V 

2r  <^(10  n=l  c12(iKn) 


dk  _  N  -1K.Z  r 

2?  RL(k)e  +  l,  A 


•v  N  K  z 

-lkz  r  .  n 


•/  _oo 


n=l 


Ln 


(2.97) 


We  have  used  eqn.  (2.14)  and  eqn.  (2.59)  in  writing  down  the  final  step. 


> 
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Equation  (2.96)  is  the  Gel' f and -Levitan  equation.  Obviously, 
it  enables  us  to  obtain  the  function  KR(x,y)  in  terms  of  the  informa¬ 
tions  on  the  reflection  and  the  transmission  coefficients. 


We  have  considered,  in  the  above  derivation,  the  situation  in 
which  a  wave  of  unit  amplitude  is  incident  from  the  left.  The  comple¬ 
mentary  problem,  in  which  a  wave  of  unit  amplitude  incident  from  the 
right  is  scattered,  can  be  dealt  with  in  a  similar  manner. 


Instead  of  eqn.  (2.77),  one  starts  with  an  analytic  function 

,  i  ■?  tv 

- - f2(x,k)e  ;  Imk  >  0  (2.98a) 


u(x,k) 


f*(X>k*)eikx 


;  Imk  <  0 


(2.98b) 


Then  one  follows  the  derivation  outlined  above. 

the  Gel ’ fand-Levitan  equation  for  K  (x,y),  viz., 

R 

roo 


ftR(x+y)  +  KR(x,y)  + 


K  (x  ,  z  )  £2  (  z+y )  dz 
K  K 

J  X 


Finally,  one  gets 


=  0  ;  y  >  x  (2 . 99) 


where 


Vz)  = 


r°°  c..(k) 

dk  11  lkz 

2tt  c  _(k)  e 
_°°  12 


-  i 


C  (iK  )  -K  Z 

11  n  n 

,  cTTTiirT  e 

n=l  12  n 


N 

I 


/•OO 


^  — OO 


(2.100) 


We  have  used  eqn.  (2.13)  and  eqn.  (2.58)  in  writing  down  the  final 
expression. 

With  this,  we  conclude  the  discussion  of  the  generalities  of 
the  inverse  problem  in  one— dimensional  quantum  scattering. 


2.3.  RECONSTRUCTION  OF  A  SYMMETRIC,  REFLECTIONLESS  POTENTIAL  IN  ONE 


DIMENSION  FOLLOWING  THE  GEL 'FAND-LEVITAN  PROCEDURE. 

Once  the  reflection  and  the  transmission  coefficients  are  known 
one  can,  in  principle,  solve  the  Gel’ f and -Levitan  equation  for  K  (x,y) 

Li 

or  K  (x,y)  and  then  obtain  V(x)  by  the  prescription  given  in  Section 
(2.2.1).  In  general,  however,  it  is  not  an  easy  task  to  get  a  solution 
of  the  Gel 1 fand-Levitan  equation  in  a  closed  form.  Only  for  a  very 
restricted  class  of  reflection  coefficients,  the  Gel ’fand-Levitan  equa¬ 
tion  is  known  to  have  simple  looking  solutions  [55,56,60].  The  simpl¬ 
est  of  the  solvable  cases  is,  perhaps,  the  one  where  the  scattering 

potential  is  reflectionless  for  all  energies.  It  is  this  simple  case 

•  « 

which  is  the  subject  of  our  investigation. 

In  what  follows,  we  will  discuss,  in  detail,  only  that  problem 
where  a  wave  of  unit  amplitude  is  incident  from  the  left.  The  results 
for  the  complementary  scattering  problem,  where  a  wave  of  unit  ampli¬ 
tude  is  incident  from  the  right,  will  only  be  quoted. 

2.3.1  Expression  for  the  Potential 

When  the  reflection  coefficient  R^(k)  has  the  simple  form 

R_^(k)  =0  V  real  k  (2.101) 

then  ft  (x,y)  in  eqn.  (2.97)  reduces  to  the  simple  form 
L 

N  K  (x+y) 

a (x,y)  =  l  Ane  n  .  (2.102) 

n=l 

Here  we  have  omitted  the  subscript  L  from  A^n  for  convenience. 

The  left  Gel ’ fand-Levitan  equation  (2.96)  then  reads 


■ 

• 

N  K  (x+y) 

l  e  +  (x ,  y  )  +  l 


n=l 


n 


N 

I 

n=l 


x 


_oo 


K^x,z)( 


<  (y+z) 


n 


dz 

(2.103) 


where  we  have  interchanged  the  summation  and  integration  on  the  right' 
hand  side.  This  integral  equation  will  accept  solutions  of  the  form 


N 


KL(x,y) 


-  I 


n=l 


g  (x)e 
n 


K  y 
n-7 


(2.104) 


where  the  functions  g  (x)  are  to  be  found.  Substituting  eqn.  (2.104) 


n 


into  eqn.  (2.103),  dividing  the  resulting  equation  throughout  by  eKn^ 
and  carrying  out  the  integration,  we  get 


N  r  K  x  N 
I  A  e  n  l 
n=l  n  v=l 


K  X 
V 


(K  +K  ) 
n  V 


K  X'i 

g^(x)  +g^(x)  +  A^e  n  J-  =  0.  (2.105) 


Hence,  the  functions  g^(x)  are  the  solutions  of  the  following  system 
of  N  linear  algebraic  equations 


K  X  N 

l 

V=1 


r  K  x 
V 

e 


n 


K  +K 

n  v 


K  X 

n 


gv(x)  +  gn(x)  +  A^e  =  0 


(2.106) 


where  n  =  l,2,...,N.  This  can  be  conveniently  written  in  the  form  of 
a  matrix  equation 


(I  +  A)g  =  C 


(2.107) 


where , 


I  =  1 1  6 


nv 


(2.108) 


A  =  A 


(k  +k  )x 
n  v 

n  (W 


(2.109) 


' 
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and 


g  = 


C  = 


r  g-j/x) 

g2(x) 


g  *  (x) 
1  N 


f  K.X'i 

-Aie 


(2.110) 


k^x 


-A2e 


K  X 
A  N 

-A  e 

v  N  ' 


(2.111) 


We  claim  that  the  coefficient  matrix  of  eqn.  (2.106)  viz.  (I+A) 
is  non-singular  and  has  a  positive  determinant  for  all  x  [49].  Proof: 
Towards  this  aim,  it  is  advantageous  to  cast  (I+A)  in  a  symmetric  form. 
Writing  eqn.  (2.106)  in  a  vector  form  and  applying  the  diagonal  trans¬ 


formation 


D  =  ||  A  1/2  6 
1 1  n  nv 


we  get  the  following  set  of  equations 


(2.112) 


(I  +  A)g  =  C 


(2.113) 


where  I  is  again  the  NxN  unit  matrix,  and 


(k  +k  )x 

_  1/2  1/2  e  n  v 

A  11  An  Av  (k  +k  ) 

n  v 


(2.114) 


g  = 


g2(x) 


(2.115) 


and 
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C  = 


1/2 

e 

1/2  K2x 
-A2.  e 


.1/2  KNX 

^  6 


(2.116) 


The  function  ^(x)  is  related  to  g^(x)  by  a  simple  relation 


£n(x)  ■  <1/2  s„(x>  • 


(2.117) 


It  can  be  immediately  checked  that 


-1 

DAD  =  A 


(2.118) 


and ,  hence 


D(I  +A)D~1  =  I  +  DAD-1  =  (I  +  A) 


(2.119) 


Obviously ,  then 


det  (I  +  A)  =  det(l+A) 


(2.120) 


So,  it  is  equivalent  to  proving  that  det(I+A)  is  positive  definite. 

A 

It  is,  for  obvious  reasons,  sufficient  to  prove  that  det  A>0 

A 

Since  A  is  a  real  symmetric  matrix,  we  know  that  it  can  always  be  dia¬ 
gonalized  [57]  i.e.  we  can  find  an  orthogonal  matrix  0  such  that 


—  1  ^  ^ 
0  AO  =  A 


(2.121) 


where  A^  stands  for  the  diagonal  counterpart  of  A.  Obviously, 

det  A,  =  det  A 
d 


So,  if  we  can  prove  that  det  A  >0,  our  initial  claim  will  be  established 

d 


It  is  well-known  that  if  a  matrix  B  is  similar  to  a  diagonal 


.. 


. 


matrix  M  then  the  diagonal  elements  of  M  are  the  eigenvalues  of  B  [57]. 
Using  this  result 


Ad  = 


X. 


0 


X. 


0 


(2.122) 


where  are  the  eigenvalues  of  A.  Clearly, 


d e t  A j  X-Xrt,*... A- T 

d  12  N 


(2.123) 


So,  det  A^  > 0  if  all  the  eigenvalues  of  A  are  positive  definite  i.e.  A 
is  a  positive  definite  matrix.  This  is  our  concern  now. 


We  know  that  a  quadratic  form  Q=Y'AY  is  positive-definite  157] 

✓\  A  ' 

if  all  the  eigenvalues  of  A  are  positive .  Hence  A  is  positive-definite  if 

(2.124) 


unless  Y  =  0 


Q  =  Y'AY  >  0 

Here  Y  is  a  N-dimensional  column  vector  and  Y’  =Y  x  .  So,  eqn.  (2.124) 
implies  that 


i.e. 


y  Y1  A  Y  >  0 
L  v  vn  n 
n,v 

(k  +k  )x 

V  -  .1/2  1/2  e  n  V 

^  n  V  (k  +K  ) 

n,v  n  v 


>  0 


i.e. 


rX 

N  1  /0  Kr,Z 

V  A 1/2  n 

)  A  ye 

L.  n  n 

J  — oo 

n=l 

dz  >  0 


(2.125) 


This  is  always  positive  unless 


N  i  to  •<  2 

1/2  -  e  "  -  0 


V  A  -*-/ 

)  A  y 

L.  n  yn 

n=l 


(2.126) 


for  all  z  between  -°°  and  x.  So,  (2.125)  is  always  positive  unless  all 


y  are  zero, 
n 


Q.  E.D. 


n ' 


. 


. 
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With  this  result,  it  is  obvious  that  the  system  of  equations 
(2.106)  has  a  unique  solution,  i.e.  g  (x)  are  unique.  The  same  is 


n 


true  for  eqns.  (2.113)  and  g  (x) . 

n 


Applying  Cramer's  rule  to  the  system  of  equation  (2.106),  we 


immediately  get  the  solution  g^(x) . 


det (I+A) 


Sn^  det  (I+A) 


R 


n  =  1,2 , . .  .  ,N 


th 


(2.127) 


where  (I+A)  is  the  matrix  with  its  n  column  replaced  by  the  vector 
R 

C  in  (2.111).  Substituting  this  value  of  g^(x)  in.  eqn.  (2.104)  we  get 


the  solution  for 


K^Cx,-) 


For  writing  down  an  expression  for  the  potential,  we  need 
K^(x,x).  From  eqn.  (2.104) 


N  K  X 

K^x,*)  =  l  gn(x)e  n 
n=l 


(2.128) 


Obviously, 


gn(x)e 


K  X 

n 


det (I+A) 
det (I+A) 


(2.129) 


where  (I+A)  is  the  matrix  (I+A)  with  its  n^  column  replaced  by  its  first 
derivative  with  respect  to  x.  Substituting  eqn.  (2.129)  in  eqn.  (2.128) 
and  recalling  the  rule  for  differentiating  a  determinant,  we  get 


K. 


_ (x,x)  =  -  - —  (£n  det(I+A)}. 
L  dx 


(2.130) 


Then  from  eqn.  (2.74),  we  have  the  potential 


V(x)  =  -2  — — 2*  Un  det  (I+A)}  . 

dx 


(2.131) 


We  can  write  eqn.  (2.131)  in  a  more  symmetric  form  using  eqn.  (2.120)  ,  viz. 


. 
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V  (x)  =  -2  -2-  {Zn  det (I+A) }  . 


(2.132) 


dx 


Had  we  started  with  R^(k)  and  the  right  Gel ' f and -Levitan  equa¬ 
tion,  we  would  have  got  [45] 


V(x)  =  -2  — 7?  {£n  det  (I+A’)} 


(2.133) 


dx 


where ,  now 


A’  = 


-(k  +k  ) X 

a1/2  1/2  e  n  v 

n  v  (k  +k  ) 
n  v 


(2.134) 


Strictly  speaking,  the  A’s  in  eqn.  (2.134)  must  have  a  subscript  R. 

The  A^n  (the  ones  in  eqn.  (2.132))  and  A^n  are  related  by  eqn.  (2.60) 
and  under  these  conditions  eqn.  (2.132)  and  eqn.  (2.133)  will  give  us 
the  same  potential  [55].  However,  in  the  case  of  a  symmetric  potential, 
which  is  of  interest  to  us,  A^n  and  ^  are  equal  and  so ,  we  prefer  to 
drop  the  subscripts  R  and  L. 


Having  got  the  potential,  we  turn  our  attention  towards  some 
of  its  properties. 

Properties  of  V(x). 

1.  V(x)  are  the  only  potentials  for  which  there  is  no  reflection. 
Proof:  We  proved  above  thdl:  the  solutions  gn(x)  °f  eqn.  (2.106)  are 
unique.  So  is,  also, K^(x, y) .  As  a  result  V(x)  is  unique  for  the  kind 
of  ^(x,y)  chosen.  And,  for  the  reflection  coefficient  satisfying  eqn. 
(2.101),  this  is  the  only  kind  possible.  This,  therefore  proves  our 
claim. 

2.  V(x)  is  negative  for  all  finite  x. 


Proof:  Let  us  write 


' 


A  =  det(A+I)  . 


(2.135) 


The  potential  V(x)  is,  then,  given  by 

v(x>  =  -  2ia-a.';.-^’.2)  . 


(2.136) 


It  is,  therefore,  sufficient  to  show  that 

V  finite  x 


aam  -  a'2  >  0 


(2.137) 


For  this  purpose,  let  us  start  by  considering  det  A.  If  we  take  out 

1/2  Kix  th  ~  1/2  Kix 

e  from  each  of  the  i  column  in  det  A  and  similarly  A.  e 

til 

from  each  of  the  i  row,  we  immediately  get 
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det  A  = 


t  N 

u  _  n 

n  a 

n=l 

e 

det  || — i — 1| 

i  n 

11  K  +K  11 

v-n=l  J 

^  n  v  J 

(2.138) 


det 


•v-s 


is  always  positive  as  it  has  the  same  form  as  det  A. 


Hence,  we  can  write 

,  ~  £x 

det  A  =  ae 

where  a  and  8  are  positive  constants.  Their  exact  appearance  is  not 

/N 

of  our  concern  now.  The  characteristic  polynomial  of  A  is  given  by 


^  N  N-l 

det{A  4-  XI}  =  X  4-  a^X  4-  ....  4- 


where  a  is  the  sum  of  the  principal  minors  of  A  of  order  r.  Putting 

A 

X  = 1  and  realizing  that  each  of  the  principal  minors  of  A  has  again  the 

A 

same  form  as  det  A,  we  immediately  get  that 


8  x 


A  =  det{A  +  l}  =  1  +  7  a  e 

u  n 
n 


(2.139) 


where  a  and  8  are  positive  constants.  Using  this  form  for  A,  we  im- 
n  n 

mediately  get 


' 
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AA"  -a’2  =  l  3^a  exp (3  x)  +y  £  a  a  (3  -3  )2exp[(3  +6  )x] 

n  n  n  l  u  n  v  n  v  n  V 

n  n,v 

(2.140) 


and  3n  being  positive,  this  gives  us  the  required  result. 

3.  V(x)  approaches  zero  exponentially  as  x  approaches  either  +°°  or  -°°. 
Proof:  Again  using  the  form  (2.139)  for  A(x)  and  (2.140)  for  AA"  -A^ 
write  down  the  potential  as 


V  (x)  = 


-2{I  exp(Bnx)  +|  I  a  av(Bn-Bv)2exp[(B  +Bu)x]} 

n  n jV 


{1+2  l  a  exp (3  x)  +  J  a  a  exp [(3  +3  )x]} 

T>  *-*  ri  \ )  *  n  \ ) 


n 


n,V 


n  V 


n  V 


(2.141) 


As  x-*-oo,  obviously,  the  numerator  goes  towards  zero  whereas 
the  denominator  tends  towards  unity.  Therefore, 

V(x)  0  ;  x  ■+  -oo  .  (2.142) 


The  behaviour  as  x-*+°°  is  implicit.  Suppose  3^  is  the  largest 

of  the  constants  3n*  The  fastest  growing  term  in  the  denominator,  as 

X++00,  is  the  one  for  which  n  =  v  =  y  i.e.  the  one  with  an  exponent 
23ux 

e  .  However,  this  term  is  absent  from  the  numerator  due  to  the  pre¬ 
sence  of  the  factor  (B  -3  )•  So,  the  denominator  grows  faster  than 

n  v 

the  numerator  as  x-»-+°°.  We  have  then 

V(x)  0  ;  x  +»  .  (2.143) 

2.3.2.  Expressions  for  the  Bound-State  Wavefunctions . 

The  simplest  way  to  get  the  expressions  for  the  normalized 
bound-state  wavefunctions  is  the  one  starting  with  the  dispersive  form 
of  u(x,k)  viz.  (2.89).  For  the  reflectionless  case,  the  last  term  on 
the  right  hand  side  of  eqn.  (2.89)  is  identically  zero.  Using  the  form 
(2.77b)  for  u(x,k)  in  the  lower  half-plane  and  using  eqn.  (2.6),  we  get 


N  f  (x,iK  )  KnX 

f  (x,-k)e  1  =1+1  -4 — r_iL_  - 

2  n=l  c12U<n)  k'1Kn 


(2.144) 


Evaluating  (2.144)  at  k = -iK  and  using  eqn.  (2.42)  and  eqn.  (2.59), 


we  get 


f2(x,iKv)  =  e 


K  X 
V 
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(k  +k  )x 
n  v 


I  A  — 7 - ; \ —  f„(x,iK  ) 

L n  (k  +K  )  2  n 

n=l  n  v 


(2.145) 


where  we  have  again  omitted  the  subscript  L  from  A^  for  reasons  al¬ 
ready  mentioned.  Multiplying  eqn.  (2.145)  through  by  A^2  and  re¬ 
membering  that  the  normalized  bound— state  wavefunction  ib  is  related 

n 

to  f  (x,iK  )  by 
2  n  J 


lpn  =  A1/2  f  (x,iK  ) 
n  n  2  n 


(2.146) 


we  have  the  following  set  of  N  algebraic  equations  for  the  N  bound- 
state  wavefunctions . 
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I 

n=l 


a1/2a1/2  1 

n  v 


(k  +k  )x 
n  v 


(k_+k  ) 
n  v 


ip  (x)  +  \p  (x) 
n  rv 


K  x 
v 

e 


(2.147) 


for  v  = 1,2 , . . .  ,N.  If  we  remember,  this  is  the  same  set  of  equations 
as  for  gn(x)  except  the  absence  of  a  negative  sign  on  the  right  hand  side. 
Applying  Cramer ' s  rule  again 


det (I+A) 

ip  (x)  =  - - —  ;  n  =  1,2,. . .  ,N  (2.148) 

n  det (I+A) 


where  (I+A)  is  the  matrix  (I+A)  with  its  n  column  replaced  by  the 
R 

column  on  the  right  hand  side  of  eqn.  (2.147).  If  we  multiply  this 
column  by  A^^e1*"1^  we  get  its  derivative.  Realizing  that 


ip  (x) 
rn 


1  det (I+A) 

A1/2eKnx  det (I+A) 
n 


n  =  1 ,2 ,  . .  .  ,N 


(2.149) 


. 


:  *2 


A  (  1^  )  A  ttl 

where  (I+A)  is  the  matrix  (I+A)  with  its  n  column  differentiated 
once  with  respect  to  x. 


Similarly,  for  the  complementary  scattering  problem,  where  a 
wave  of  unit  amplitude  is  incident  from  the  right,  the  expression  for 
^n(x)  is  145] 


VX> 


1 


A1/2 


det(H-A’) 
det (I+A* ) 


(2.150) 


^  ~  ('n') 

where  the  matrix  A'  was  introduced  in  eqn.  (2.134).  (I+A*)  carries 

A  (n) 

the  same  meaning  as  (I+A)  .  Again,  we  have  dropped  the  subscript 

R  from  A  in  eqn.  (2.150)  for  reasons  already  discussed. 


2.3.3.  Imposing  Symmetry  on  the  Potential. 

We  have  been  able  to  derive,  so  far,  the  expressions  for  cal¬ 
culating  a  transparent  potential  V(x)  and  the  corresponding  bound- 
state  wave  functions  ipn(x)  i-n  ternis  of  2N  parameters,  N  being  the 
number  of  bound-states.  The  N  k_^ts  are  the  bound-state  energies.  The 
other  N  parameters,  viz.,  the  A_^'s,  have  to  be  fixed  by  some  explicit 
information  regarding  the  bound-state  wavefunctions  or  some  desired 
constraint  on  the  potential  V(x).  We  are  specifically  interested  in 
a  potential  which  is  symmetric  about  the  origin,  i.e. 

V (x)  =  V(-x)  .  (2.151) 


We  claim  that  the  condition  (2.151)  is  guaranteed  provided  that  the 
A^'s  are  given  by 


n 

m^n 


K  +K 

m  n 


K  -K 

m  n 


m  =  1,2, _ ,N  (2.152) 


Proof:  We  have  already  seen  that 


. 


(2.153) 


det{ A  + i)  =  1  +  +  ....  +  a^  ^  +  a^ 


where  a^_  is  the  sum  of  the  principal  minors  of  A  of  order  r.  Also, 

A  A 

det  A  is  given  by  eqn.  (2.138).  All  the  principal  minors  of  A  have 
the  same  form  as  (2.138)  except  that  one  or  more  (n,v)  are  missing 
depending  upon  the  order  r. 

N 

n5lKnx 

Let  us  take  out  a  factor  e  on  the  right  hand  side  of  eqn. 

(2.153)  yielding 

N  N  N 

y  k  x 

-  n=l  n 

det(A  +  l)  =  e 


- 1 
n=l 
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+  e 
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N 
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K  X 
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(2.154) 


Here,  for  example, 


and 


N  A  2k 


li-trA-  l 


i  2K 
P  =  1  p 


(2.155) 


aN-l  =  d6t 


(  1 
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r  N 

II  1  1 
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i  i  n 
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11  K  +K  1 
^  n  v 
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(n,v)=p 
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2  y  k  x 

SPp  n 


missing 


n^p 


(2.156) 


a.T  =  det  A  . 
N 


(2.157) 


From  eqn.  (2.132)  it  is  clear  that  the  overall  exponential  in 
eqn.  (2.154)  does  not  contribute  to  the  potential.  A  symmetric  poten¬ 
tial  then  requires  that  the  quantity  in  the  braces  of  eqn.  (2.154) 
be  symmetric.  It  is  important  to  note  that  the  only  dependence  on  x 


•  9 

. 

is  in  the  exponents.  So,  the  coefficients  of  exponents  differing  only 
in  an  overall  sign  have  to  be  the  same. 


Let  us  consider  the  first  term  in  the  braces  in  eqn.  (2.154) 


The  exponent  is 
N 

-  L  Knx 
n=l 


with  a  factor  of  unity.  So,  unity  must  also  be  the  coefficient  of 
N 

1  Knx 
„n=l  n 


for  symmetry.  This  comes  from  a  . 
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n  V 


all 

n,v 


f  N 

II.  A 

n=l  n 


We  must  have,  then 


=  1 


(2.158) 


According  to  Lemma  1  in  Appendix  A,  the  condition  (2.158)  becomes. 
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(2.159) 


Similarly,  the  p  =  i  term  in  a,  gives  us  the 


coef f . of  e 


(_K1  "K2  ” - "Ki-1+Ki  "Ki+1  " - 'KN)x 
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Its  symmetric  counterpart  in  a^_^  gives  us ,  the 
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missing 


(2.160) 


(2.161) 


Again  using  the  Lemma  mentioned  above,  the  equality  of  these  two 


coefficients  dictates  to  us 
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(2.162) 

where  IT’  implies  that  the  factor  with  j=i  is  to  be  excluded  in  the  product. 


Multiplying  both  the  sides  of  eqn.  (2.162)  by 
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and  using  condition  (2.159)  on  the  right-hand  side  of  the  resulting 


equation,  we  get 
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(2.164) 


Remembering  that  A^'s  have  to  be  positive,  eqn.  (2.164)  gives 

;  j ,i  =  1,2, . . . ,N 


A. 

i 
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which  is  the  desired  result. 


Let  us  now  try  to  prove  the  converse,  i.e.,  if  eqn.  (2.152) 
holds  good  then  V(x)  is  symmetric  [47].  For  this  purpose,  let  us  cast 

eqn.  (2.153)  in  a  slightly  different  form.  We  observed  earlier  that 

^  /\ 

all  the  principal  minors  of  det  A  have  the  same  form  as  det  A  except 
that  some  of  the  rows  and  columns  are  missing.  If  we  define  by  s,  all 
the  subsets  of  {l,2,...,N}  including  the  null  set  and  the  full  set,  we 


■ 


' 


can  rewrite  eqn.  (2.153)  as 
det{A  +  I>  =  l  a 


(2.165) 


A 

where  a  is  a  principal  monor  of  A  with  only  those  rows  and  columns 

O 

which  belong  to  s.  We  realize  that  a^  corresponding  to  the  null  set 
is  just  unity.  Obviously,  then 
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(2.166) 


n,m£s 


Using  Lemma  1  of  Appendix  A  again,  the  determinant  on  the  right  hand 
side  of  eqn.  (2.166)  can  be  written  as 


K  X 


a  = 
s 


( 

l  ad 

r 

e 

p£S 

p  ' 

^p£s 

V 

* 

N 

o 

N 

n  (Kn-K 

,)2 

n 

(K  - 

i>l  1 

J 

z 

n 

^p£s 


2k 


N 


i  >N-1 


N-l  j 


N  2  N  2  N  2 

n  (K..-HO  n  (k  +k.)  _  n  (k  ,+k.) 

j>l  1  J  j>2  2  2  j  >N-1  1,-1  2 


.  (2.167) 

all  K  ' s  £s 
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The  factor  in  brackets  contains  only  those  k^’s  i-n  t^ie  product  which 
are  in  S.  Others  are  missing.  Supplementing  this  with  Lemma  2  in 
Appendix  A,  we  can  write  it  in  a  more  intelligible  form 
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(2.168) 


Substituting  eqn.  (2.168)  in  eqn.  (2.165),  we  have 
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(2.169) 


If  we  now  substitute  the  assumed  conditions  for  (A  /2k  ) ,  we  get 

P  P 
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(2.170) 


Obviously,  those  terns  in  the  first  parentheses  which  have  n£s  cancel 
out  with  the  ones  in  the  second  parentheses.  We  are,  thus,  left  with 

2  y  k  x 

^  u  p 

det{A  +  l}  =  l  e  p£S  II(s,s)  (2.171) 
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(2.172) 


and  s  stands  for  the  complement  of  s.  It  is  clear  that  II(s,s)  is  sym¬ 
metric  under  the  interchange  of  s  and  s,  i.e., 

n(s,s)  =  n(I,s).  (2.173) 


We  can  write  eqn.  (2.171)  as: 
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(2.174) 


Since  II(s,s)  is  symmetric  under  the  interchange  s  + s  and  the  sum  runs 
over  all  possible  subsets,  we  can  write 
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.  (2.175) 


The  first  term,  as  before,  does  not  affect  the  potential  and  the  second 


term  is  manifestly  symmetric  in  x.  So,  the  potential  is  symmetric.  QED. 


2.3.4.  A  Numerical  Test  of  the  Technique. 

As  a  numerical  test  of  the  above  technique,  let  us  consider  the 
Poschl-Teller  potential 

f 

V(x)  =  -  -(-A~P  (2.176) 

cosh  x 

r 

with  A>1.  This  is,  obviously,  asymmetric  potential.  This  potential 
has  been  studied  and  is  known  to  be  an  exquisite  example  of  a  reflec¬ 
tionless  potential  for  certain  values  of  A  [61].  The  bound-state 
energies  of  (2.176)  are 

E  =  -(A  -1  -n)2  ;  n  <  A  -  1  (2.177) 

n  — 

and  n  can  take  the  values  0,1,2,....  Also  the  reflection  coefficient 
is  known  to  be 

|R^(k)|2  =  cos^(cj)^  -  ({O  (2.178) 


where 
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Obviously,  for  all  integer  values  of  A (=1,2, 3, . . .)  the  reflection  co¬ 
efficient  identically  vanishes  for  all  real  k. 

So,  this  is  obviously  the  perfect  case  for  a  test  of  our  method. 
Our  technique  should  reconstruct  potential  (2.176)  for  integer  values 
of  A  if  we  use  the  bound-state  energies  given  by  eqn.  (2.177).  And, 
in  fact,  it  does.  We  reconstructed  the  Poschl-Teller  potentials  with 
A  =  2, 3, 4, 5  and  6.  It  was  exact  in  each  case.  We  display  the  result 
for  A  =  6  in  Figure  2.6. 

That  proves  the  uniqueness  and  credibility  of  the  technique 


outlined  in  this  chapter. 


VRLUE 
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Fig.  2.6.  Results  for  the  Poschl-Teller  potential  (A -6): 

(a)  Actual  (A)  and  the  reflectionless  (E)  potentials:  exact  agreement. 

(b)  Bound-state  wavefunct ions . 


3.  APPLICATIONS  TO  SOME  SYMMETRIC  POTENTIALS  IN  ONE  DIMENSION: 

RELIABILITY  OF  THE  APPROXIMATION. 

Given  a  sequence  of  bound-state  energies,  we  can  uniquely  con¬ 
struct  a  symmetric,  reflectionless  potential  which  supports  those 
bound-states.  This  technique  does  not  have  much  appeal  as  an  exact 
inverse  method.  There  are  not  very  many  reflectionless  potentials  of 
interest.  Moreover,  we  do  not  know,  a  priori,  whether  a  given  sequence 
of  bound-state  energies  corresponds  to  a  reflectionless  potential  or 
not.  What  we  are  interested  in  is  the  reliability  of  this  technique 
as  a  method  for  approximately  constructing  a  symmetric  potential  from 
its  bound-state  energies.  In  other  words,  if  we  are  given  the  bound- 
state  energies  of  a  symmetric  potential,  which  is  not  necessarily  re¬ 
flectionless,  then  how  closely  the  constructed  reflectionless  potential 
resembles  the  actual  one.  The  obvious  way  to  check  this  is  to  apply 
this  technique  to  some  known  symmetric  potentials  in  one  dimension. 

This  is  the  aim  of  this  chapter. 

Section  3.1  is  devoted  exclusively  to  non-confining  potential 
wells,  Section  3.2  to  confining  potentials  and  in  Section  3.3  we  dis¬ 
cuss  the  problem  of  convergence. 

mt 

3 . 1  Non-Confining  Potential  Wells 

In  this  class,  we  consider  only  those  symmetric  potential  wells 
which  are  negative  and  non-singular  on  the  infinite  interval  (-°°<x<00) 
and  which  decay  fast  enough,  as  x  =  ±°°,  to  support  only  a  finite  number 
of  bound  states.  In  all  the  cases,  we  use  the  exact  bound-state  ener¬ 
gies  to  construct  the  corresponding  reflectionless  potential  and  the 
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estimated  bound-state  wavefunctions .  We  also  calculate  the  reflection 
coefficient  and  plot  it  as  a  function  of  the  momentum  ’k’. 

3.1.1.  Poschl-Teller  Potentials  with  Non-Integer  A's. 

The  Poschl-Teller  potential  given  by  eqn.  (2.176)  is  not  re¬ 
flectionless  for  a  non- integer  value  of  A.  This  potential  satisfies 
all  the  properties  listed  above  and  its  reflection  coefficient  displays 
a  very  simple  behaviour  with  increasing  momentum  k. 

Using  the  exact  bound-state  energies,  reflectionless  potentials 
corresponding  to  Poschl-Teller  potentials  with  A  =  (2.1  to  2.9)  and 
A  =  (3.1  to  3.9)  were  calculated.  In  Figure  3.1  to  Figure  3.5,  we 
have  summarized  the  calculations  for  A  =  3.1  to  3.5.  These  five  cases 
suffice  to  elucidate  the  features  of  the  reflectionless  approximation. 

As  we  see  from  eqn.  (2.178)  and  eqn.  (2.179),  the  reflection 
coefficient  for  a  non-integer  A  decays  monotonically  with  increasing 
k.  At  a  value  A  =  (n  +  0.5);  n  =  l,2,...,  the  rate  of  decay  is  the 
slowest.  The  reflection  coefficient  being  periodic  in  A,  as  A  varying 
both  ways  from  (n  +  0.5)  to  an  integer  value  i.e.  towards  n  or  (n  + 1)  , 
the  reflection  coefficient  shows  exactly  the  same  behaviour.  The  closer 
A  is  to  the  integer  value,  the  faster  the  reflection  coefficient  decays. 
And,  as  is  obvious  from  the  results,  the  faster  the  reflection  coeffi¬ 
cient  decays  with  increasing  k  the  better  is  the  agreement  between  the 
reflectionless  approximation  and  the  actual  potential. 

To  further  elucidate  this  behaviour  the  Poschl-Teller  potential 
with  A  =  3.0001  was  reconstructed.  As  we  see  in  Figure  3 . 6  ,  the  reflec¬ 
tion  coefficient  is  equal  to  unity  at  k = 0  and  practically  zero  every- 
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(a)  Actual  (A)  and  the  reflectionless  (E)  potential. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Results  for  the  Poschl-Teller  Potential  (X  =  3.3): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Results  for  the  Poschl-Teller  Potential  (A  =  3.5): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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where  else.  The  reflectionless  potential,  though  not  the  same  exactly, 
is  extremely  close  to  the  actual  potential. 

The  potentials  with  X  =  2.5,  3.5,  4.5  and  5.5  have  exactly  the 
same  reflection  coefficient.  It  was,  therefore,  interesting  to  see  if 
there  is  a  difference  in  the  agreement  between  the  reflectionless  and 
the  actual  potential  in  these  cases.  We  see  from  Figure  3. 7  and  Figure 
3.8  that  the  overall  agreement  gets  better  as  the  potential  gets  deeper 
despite  the  fact  that  the  reflection  coefficient  remains  the  same. 

3.1.2.  The  Gaussian  Well. 

Second  in  this  category  was  a  well  of  the  kind 

2 

V (x)  =  — V  eX  ,  _oo  <  x  <  oo  (3.1) 

For  obvious  reasons,  we  call  it  a  Gaussian  Well.  V  is  the  depth  of 

o 

the  well.  Five  such  wells  with  different  V  were  considered. 

o 

Since  the  analytic  solution  to  the  Schrodinger  equation  with 
such  a  potential  is  not  known,  the  eigenvalues  and  the  reflection  co¬ 
efficient  in  each  case  were  calculated  using  the  numerical  methods  of 
Appendix  B.  Eigenvalues  and  the  zeros  of  the  reflection  coefficients 
have  been  tabulated  in  Table  3.1  and  Table  3.2.  Reflection  coefficients 
and  the  results  of  the  inverse  calculations  are  displayed  in  Figure  3.9 
to  Figure  3.13. 

A  zero  of  the  reflection  coefficient  is  that  value  of  momentum 
where  the  reflection  coefficient  is  zero  to  five  places.  Since  we  were 
interested  in  the  general  systematics  of  the  reflection  coefficient. 
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Table  3.1  Eigenvalues  of  the  Gaussian  Well 

-x2 


V (x)  =  -V  e 
o 


—o°  <  x  <  +°° 


Ns\vEig  envalues 
xin  ascend- 
\ing 

Depth  \order 

V  \ 

o  \ 

_  2 

E1  =  -K1 

E  =  -K2 

2  2 

E3  "  -K3 

E4  =  -K4 

E  =  -K2 
b5  5 

2.0 

(-0.955)* 

(-0.954)f 

• 

« 

00 

• 

o 

(-5.549)* 

(-5.549)+ 

(-1.568)* 

(-1.568)+ 

18.0 

(-14.135)* 

(-14.135)+ 

(-7.256)* 
(-7 .256) ^ 

(-2.278)* 

(-2.278)1* 

30.00 

(-24.900)* 

(-24.900)+ 

(-15.525)* 

(-15.525)+ 

(-7.918)* 

(-7.919)+ 

(-2.429)* 

(-2.430)+ 

46.0 

(-39.595)* 

(-39.594)+ 

(-27.592)* 

(-27.592)+ 

(-17.289)* 

(-17.290)+ 

(-8.901)* 

(-8.903)+ 

(-2.809)* 

(-2.810)+ 

Note:  *  Corresponds  to  the  numerical  solution  of  the  eigenvalue  equation, 

t  Corresponds  to  the  matrix  method. 

These  symbols  carry  the  same  meaning  in  the  following  tables  of 
eigenvalues,  too. 


Table  3.2 


Zeros  of  the  Reflection  Coefficient:  The  Gaussian  Well 


2 

V(x)  =  -V  e  X  ;  -oo  <  x  <  oo 
o 


Momentum  at 
^\which  R  is 

Depth 

V 

o 

ist 

2nd 

2.0 

1.04 

3.46 

00 

• 

o 

0.82 

3.52 

18.0 

3.54 

30.0 

0.28 

3.80 

46.0 

3.92 

VALUE 
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Fig.  3.9.  Results  for  the  Gaussian  Well  (VQ  =  2.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.10.  Results  for  the  Gaussian  Well  (VQ  =  8.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.11.  Results  for  the  Gaussian  Well  (VQ  =  18.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.12.  Results  for  the  Gaussian  Well  (VQ  =  30.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.13.  Results  for  the  Guassian  Well  (VQ  =  46.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


no  attempt  was  made  to  obtain  the  zeros  to  a  higher  accuracy.  As  we 
shall  see  shortly,  these  values  prove  sufficient  for  our  purpose. 

The  reflection  coefficient,  in  this  case,  has  a  complex  be¬ 
haviour.  In  the  absence  of  an  analytic  form  for  it,  it  is  difficult 
to  predict  anything  about  its  behaviour  for  a  general  'kT.  However, 
in  the  range  of  momentum  considered,  viz.  k  =  0.0  to  5.0,  it  shows  some 
regular it ies. 

Firstly,  it  decays  very  fast  and  has  a  few  zeros.  Secondly,  it 
shows  certain  periodicity  with  increase  in  V  .  For  instance,  the  be¬ 
haviour  of  the  reflection  coefficient  for  V  =2.0,  .8.0  and  30.0  are 
similar  -  in  each  case,  it  has  a  zero  close  to  k  =  0.  However,  for 
Vo  =  18.0  and  46.0,  it  has  a  zero  at  a  very  large  k  where  its  amplitude 
is  practically  zero  anyway. 

The  agreement  between  the  reflectionless  approximation  and  the 
actual  potential  is,  once  again,  very  good.  The  fast  decay  of  the  re¬ 
flection  coefficient  with  increase  in  k  again  seems  to  be  a  favourable 
factor. 


Within  the  scope  of  these  numerical  results,  it  seems  difficult 
to  go  beyond  these  naive  observations. 

3.1.3.  The  Harmonic  Oscillator  Well. 

As  the  next  application,  the  following  potential  was  considered 


fo  ;  |x |  >  a 

-E  +x^  ;  1  x  1  <  a 

V  O  ~ 

(3.2) 

Obviously, 


•  ' 


74 


E  =  a2  .  (3.3) 

o 

An  increase  in  E^ ,  thus,  increases  both  the  depth  and  the  range  of  the 

potential.  Five  such  potentials  with  different  E  were  considered. 

o 

The  eigenvalue  equations  for  such  potentials  have  already  been 
derived  in  Appendix  B.  In  order  to  calculate  the  eigenvalues  of  this 
potential  using  eqn.  (B.10)  and  eqn.  (B.ll)  we  just  need  to  know  the 
even  and  odd  solutions  of  the  Schrodinger  equation  in  the  region 
-a<x<_a.  The  Schrodinger  equation  in  this  region  looks  like 

2 

+  [E  +  a2  -  x2]ij;  =  0  .  (3.4) 

dx 

If  we  change  the  variable  to 
y  =  Jl  x 

this  equation  is  transformed  into 

-  (t-  y2  +  a)ip  =  0  (3.6) 

dy 


where 

a  =  -  \  (E  +  a2) 


(3.7) 


Equation  (3.6)  is  the  parabolic  cylinder  equation  [62].  The  even  and 
odd  solutions  are  given  by  [62] 


ip  (x,k)  =  A  e 
e  1 


-W 


A  ,  1 

M(y  a  + 


1 

2’  X  } 


(3.8) 


and 


\ \>  (x,k)  =  A_xe 
o  z 


-^x  A  ,3  2  2 . 

M(j  a  +  4,  J »  x  ) 


where  A^  and  A^  are  the  normalization  constants  and  M(a,b,x) 


(3.9) 
is  the 


•  -  ' 
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confluent  hyper geometric  function. 

These  solutions,  for  E  <0,  were  substituted  in  eqn.  (B.10)  and 
eqn.  (B.ll)  of  Appendix  B  to  obtain  the  eigenvalue  equations  for  this 
Particular  problem.  The  resulting  equations  were,  then,  numerically 
solved  by  the  Regula-Falsi  technique  to  obtain  the  bound-state  energies 
[63]. 


Similarly,  for  E  >0,  these  solutions  were  used  in  eqn.  (B.14) 
of  Appendix  B  to  calculate  the  reflection  coefficient  for  various  ener¬ 
gies. 

In  both  of  these  calculations ,  values  of  some  confluent  hyper- 
goemetric  functions  M(a,b,x)  were  required  for  some  values  of  the  para¬ 
meters  a,  b  and  x.  They  were  calculated  by  summing  their  respective 
infinite  series  [62]  to  a  sufficient  number  of  terms. 

Table  3.3  and  Table  3.4  list  the  eigenvalues  and  zeros  of  the 
reflection  coefficient,  respecitvely ,  for  each  of  the  five  potentials 
of  this  kind  which  were  considered.  The  inverse  calculations  and  the 
reflection  coefficients  are  plotted  in  Figure  3.14  to  Figure  3.18. 

The  behaviour  of  the  reflection  coefficient  as  a  function  of  k 
is  reminiscent  of  the  diffraction  pattern  in  optics.  Borrowing  the 
terminology  from  there,  the  reflection  coefficient  shows  a  principal 
maximum  and  some  pronounced  secondary  maxima.  The  agreement  between 
the  reflectionless  and  the  actual  potential,  though  tolerable,  is  not 
as  good  as  in  the  case  of  Poschl-Teller  potentials  or  the  Gaussian  wells. 


The  secondary  maxima  and  the  overall  amplitude  of  the  reflection 


* 


. 


•  rj 


■ 


Table  3.3  Eigenvalues  of  the  Harmonic  Oscillator  Well 


V  (x) 


>  a 


with 


<  a 


E 

o 


2 


=  a 


\Eigenvalues 
\  in  Ascen- 
N.  ding 

E  \order 

o 

Ei  ■ 

E2  *  -K2 

E  =-K2 

3  K3 

E4  '  -K4 

E  =  -K2 

E5  5 

2.0 

-1.069 

(-1.069)* 

(-1.069)+ 

4.0 

-3.006 

(-3.006)* 

(-3.005)+ 

-1.092 

(-1.092)* 

(-1.091)+ 

6.0 

-5.001 

(-5.001)* 

(-5.000)+ 

-3.011 

(-3.011)* 

(-3.011)+ 

-1.106 

(-1.106)* 

(-1.106)+ 

8.0 

-7.000 

(-7.000)* 

(-7.000)f 

-5.002 

(-5.002)* 

(-5.001)+ 

-3.016 

(-3.016)* 

(-3.016)+ 

-1.117 

(-1.117)* 

(-1.117)+ 

10.0 

-9.000 

(-9.000)* 

(-9.000)+ 

-7.000 

(-7.000)* 

(-7.000)+ 

-5.002 

(-5.002)* 

(-5.003)+ 

-3.020 

(-3.020)* 

(-3.021)+ 

-1.125 

(-1.125)* 

(-1.126)+ 

Note:  The  numbers  without  parentheses  correspond  to  the  calculation 

where  the  confluent  hypergeometr ic  functions  were  used  in  the 
eigenvalue  equations. 
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Table  3.4  Zeros  of  the  Reflection  Coefficient:  The  Harmonic  Oscillator 

Well. 


V  (x) 


fo 

1  x  1  >  a 

> 

—  with  E 

-E  +x^  ; 

^  o 

o 

v  I 

X 

2 

a 


Momentum 
\at  which 
R  is 

E  \zero 
o  \ 

ist 

2nd 

3rd 

4th 

5  th 

6th 

2.0 

0.7405 

2.425 

3.65 

4.85 

4.0 

0.7027 

2.12 

3.10 

3.95 

4.80 

6.0 

0.6854 

1.9992 

2.867 

3.65 

4.35 

8.0 

0.6748 

1.9318 

2.75 

3.45 

4.15 

4.75 

10.00 

0.6673 

1.8876 

2.6638 

3.329 

3.95 

4.55 

* 


VALUE 

2.0  2-4  2.0  3.2  3.6  4.0 
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Fig.  3.14.  Results  for  the  Harmonic  Oscillator  Well  (EQ  =  2.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 


X  VALUE 

0.0  0.4  0.0  1.2  t.S  2.0  2.4  2.0  3.2  3.8  4.0 
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Fig.  3.15.  Results  for  the  Harmonic  Oscillator  Well  (E0  =  4.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


X  VALUE 

0.0  0.4  0.0  1.2  1.1  2.0  2.4  2.0  3.2  3.0  4.0 
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Fig.  3.16.  Results  for  the  Harmonic  Oscillator  Well  (EQ  =  6.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


X  VALUE 

0.0  0.4  0.0  1.2  l.s  2.0  2.4  2.0  3.2  3.0 
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Fig.  3.17.  Results  for  the  Harmonic  Oscillator  Well  (EQ  =  8.0): 

(a)  Actual  (A)  and  the  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  for  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


X  VALUE 

0.0  0.4  0.0  1.2  1  .6  2.0  2.4  2.6  S.2  3.6  4.0 
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Fig.  3.18.  Results  for  the  Harmonic  Oscillator  Well  (E0  =  10.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


coefficient  get  higher  as  E  increases.  The  successive  zeros,  however, 
seem  to  be  getting  closer  to  each  other.  The  change  in  the  first  zero 
is  nominal  but  the  others  seem  to  cluster  faster. 

The  inverse  reconstructions  show  that  the  reflectionless  poten¬ 
tials  tend  to  settle  down  to  a  stable  pattern,  which  is  not  very  close 

to  the  actual  potential,  with  increase  in  E  .  This  observation  is,  of 

o  * 

course,  limited  to  the  five  cases  which  we  have  considered.  There  is, 
however,  no  suggestion  either  that  the  reflectionless  potential  will 
eventually  diverge  from  the  actual  potential.  So,  we  will  still  clas¬ 
sify  this  case  as  a  favourable  one. 


3.1.4.  The  Linear  Well. 

We  consider  the  potential  well  given  by 


V(x) 


fO 

-E  +  lx 
o  1 


Here 


E 

o 


=  a 


x|  _>  a 
x|  <_  a 


(3.10) 


(3.11) 


Again,  a  change  in  Eq  changes  both  the  depth  and  the  range  of  the  po¬ 
tential  well.  Five  such  potentials  with  different  E  were  considered. 

o 

Our  first  concern,  as  usual,  is  to  calculate  the  exact  eigen¬ 
values  and  reflection  coefficient  of  each  of  these  potential  wells.  In 
this  case,  the  numerical  methods  of  Appendix  B  were  used.  Table  3.5 
and  Table  3.6  summarize  the  results  of  these  calculations. 


Inverse  calculations  and  the  reflection  coefficients  have  been 


plotted  in  Figure  3.19  to  Figure  3.23. 


. 
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Table  3.5  Eigenvalues  of  the  Linear  Well 


1° 

a 

l-E0+|x| 

1  V  | 

X 

•  0S 

with  E  =  a 
a  o 

N.  Eigenvalues 
in  Ascen- 
\ding  or- 

E  \der 

o 

E  =-K2 
hl  K1 

E2  =  -K2 

E  =-K2 
h3  K3 

E4=-K4 

E  -  -K2 

5  5 

1.67845019 

(-0.692)* 
(-0.691) f 

2.79315249 

(-1.775)* 

(-1.775)+ 

(-0.491)* 

(-0.490)+ 

3.66807351 

(-2.649)* 

(-2.649)+ 

(-1.332)* 

(-1.332)+ 

(-0.454)* 

(-0.454)+ 

4.45402433 

(-3.435)* 

(-3.435)+ 

(-2.116)* 

(-2.116)+ 

(-1.209)* 

(-1.209)+ 

(-0.402)* 

(-0.402)+ 

5.17032952 

(-4.152)* 
(-4. 151) + 

(-2.832)* 

(-2.832)+ 

(-1.922)* 

(-1.922)1' 

(1.086)* 

(1.086)+ 

(-0.386)* 

(-0.385)f 

Table  3.6  Zeros  of  the  Reflection  Coefficient:  The  Linear  Well 
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VALUE 

3,0  3.6  4 .2  4.6  6;4  6.0 
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Fig.  3.19.  Results  for  the  Linear  Well  (EQ  =  1.67845019): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


X  VALUE 

0.0  0.6  !.?  1  .1  2.4  3.0  3.6  4.2 
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Fig.  3.20.  Results  for  the  Linear  Well  (EQ  =  2.79315249): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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X  VALUE 
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Fig.  3.22,  Results  for  the  Linear  Well  (E0  =  4.45402433): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.23.  Results  for  the  Linear  Well  (E0  =  5.17032952): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 
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Despite  its  overall  resemblance  to  the  reflection  coefficient 
of  the  harmonic  oscillator  well,  the  reflection  coefficient  of  the 
linear  well  has  some  very  distinctive  features. 

The  number  of  observable  peaks  (up  to  five  decimal  places) 

shows  a  fast  proliferation  with  increase  in  E  in  contrast  to  the  har- 

o 

monic  oscillator  well.  However,  the  rate  of  decay  of  the  reflection 
coefficient  as  a  function  of  momentum  is  much  faster  in  this  case  as 
compared  to  the  harmonic  oscillator  well.  With  these  two  features 
together,  as  one  would  expect,  the  peaks  are  narrower  in  the  case  of 
the  linear  well. 

For  each  E  ,  the  gaps  between  the  successive  zeros,  in  the  case 
of  the  harmonic  oscillator  well,  steadily  decreased.  In  this  case, 
they  show  some  kind  of  periodicity.  The  first  gap  is  larger,  the 
second  is  smaller  and  so  on. 

The  corresponding  gaps  in  the  zeros  for  two  different  Eq’s  are 
again  periodic.  With  Eq =  1. 67845019 ,  for  instance,  the  first  is  smaller, 
the  second  is  larger  and  so  on.  For  E^  = 2 . 79315249 ,  on  the  other  hand, 
the  first  is  larger,  the  second  is  smaller  and  so  on.  This  cycle  alter¬ 
nates  with  each  E  . 

o 

The  peaks  do  get  narrower  with  increasing  E^  but  not  necessarily 
in  the  same  order  because  of  the  alternating  behaviour  described  above. 

Now,  if  we  look  at  the  reflectionless  approximation  we  can  im¬ 
mediately  discern  a  better  agreement  in  this  case  between  the  reflec— 
tionless  and  the  actual  potential  as  compared  to  the  agreement  in  the 
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case  of  the  harmonic  oscillator  well.  First  obvious  reason  seems  to 
be  the  fact  that  the  reflection  coefficient  decays  much  faster  in  this 
case. 

Narrow  peaks  seem  to  favour  the  agreement .  This  was  also  ob¬ 
vious  in  the  case  of  the  harmonic  oscillator  well  in  a  benign  way. 

It  is  much  more  pronounced  here  as  the  overall  decay  does  not  seem  to 
slow  down  appreciably  with  increase  in  . 

The  periodicity  in  the  behaviour  of  the  reflection  coefficient 
with  increasing  seems  to  be  reflected  in  the  periodicity  of  the 
agreement  between  the  two  potentials.  The  dependence  is  curious  but 
it  cannot  be  unravelled  by  such  a  numerical  exercise. 

This  potential  seems  to  be  much  more  promising  than  the  harmonic 
oscillator  well  from  the  point  of  view  of  this  approximation. 


3.1.5.  The  Finite  Square  Well. 

We  now  come  back  to  the  most  familiar  potential  of  elementary 
quantum  mechanics,  viz.,,  the  finite  square  well 


0  ;  x  <  -a 


V  (x) 


-a  <  x  <  a 
x  >  a 


(3.12) 


For  E<0,  we  get  bound  states.  Defining 

k2  =  (E  +  E  )  =  (E  -  |  E  | ) 
o  o 

K2  =  -E  =  I E  I 


(3.13) 

(3.14) 


the  eigenvalue  equation  for  the  even— parity  levels  is  given  by  [61] 


. 


■ 
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k  tan  ka  =  K 

and  for  the  odd-parity  levels  is 
k  cot  ka  =  -K 

The  even  eigenfunctions  are  given  by 


A  cos  k  x 
+  n 


ip  (x)  =  - 
en 


;  0  <  x  <.  a 


K  (a-x) 

,  .  n 

A,  cos  k  a  e  ;  x  >  a 

+  n  ’ 


ip  (-x)  =  ip  (+x) 
en  en 


(3.15) 


(3.16) 


(3.17) 


Similarly,  the  odd  eigenfunctions  are  given  by 

;  0  <  x  <  a 


\p  (x)  =  ■ 
on 


A  sin  k  x 
-  n 


K  (a-x) 

.  .  .  n 

A  sm  k  a  e  ;  x  >  a 

-  n  3 


ip  (-x)  =  -  ip  (x) 
^  on  on 


The  normalization  constants  are  given  by 


A+  n 


for  both  cases. 


(3.18) 


(3.19) 


For  E>0,  we  get  the  continuum  states.  The  reflection  for  this 
case  is  well-known  [64] 


R  = 


1  + 


4 1 E 

( 

111 

\ 

+  E0) 

Eq  Sin2{2 J 

(W 

+  eoa 

A  JJ 

-1 


(3.20) 


where 


a 


(3.21) 


. 


Five  such  potential  wells  were  considered.  All  had  the  same 

7T 

width,  viz.,  a=  —  .  They  were,  however,  of  different  depths.  Eigen¬ 
values  were  calculated  by  solving  eqn.  (3.15)  and  eqn.  (3.16)  numeri¬ 
cally  by  the  Newton-Raphson  method  [63].  Reflection  coefficients  were 
found  by  using  eqn.  (3.20).  Table  3.7  and  Table  3.8  list  the  eigen¬ 
values  and  the  zeros  of  the  reflection  coefficients  respectively. 

jm 

The  results  of  the  inverse  calculations  and  the  reflection  co¬ 
efficients  are  plotted  in  Figure  3.24  to  Figure  3.28. 


This  is  a  highly  reflecting  potential.  For  smaller  depths, 
the  reflection  coefficient  is  relatively  smaller,  but  still  higher  than 
the  harmonic  oscillator  well,  or  the  linear  well.  With  an  increase  in 
E^,  the  peaks  get  wider  and  higher.  In  fact,  the  reflectionless  ap¬ 
proximation,  though  tolerable  for  wells  of  smaller  depths,  disagrees 
considerably  with  the  actual  potential  as  E^  is  increased.  There  is  no 
indication,  whatsoever,  that  it  converges  to  the  actual  potential.  It 
does  not  even  seem  to  settle  down  to  a  stable  pattern;  thus,  exhibiting 
a  diverging  trend.  We,  therefore,  classify  this  potential  among  the 
unfavourable  ones  for  our  approximation. 


3.1.6.  The  Secant-Square  Well. 

Finally,  we  consider  a  potential  well  of  the  type 


0 


V(x) 


8 


cos^2x 


x  |  _>  a 
x|  <  a 


(3.22) 


Obviously  the  range  a  and  E^  are  related  by 
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Table  3.7  Eigenvalues  of  the  Finite  Square  Well 


V  (x) 


0  ;  x  <  -  tt  /  2 

I  -E0  ;  -tt/2  <  x  <  tt/2 

[o  ;  x  >  tt/2 


NEig enva lu  e  s 
\  in  Ascen- 
\  ding 

E  \rd" 
o  \er 

E  —  £2 
*•1  _K1 

E  =  -K2 

fc2  K2 

E  =  -K2 

3  3 

„  2 
E4  =  '"A 

„  2 
E5  =  -K5 

E6  =  ~Kl 

2.5 

-2.002 

(-2.002)* 

-0.680 

(-0.680)* 

6.5 

-5.864 

(-5.864)* 

-4.017 

(-4.017)* 

-1.257 

(-1.257)* 

12.5 

-11.784 

(-11.784)* 

-9.664 

(-9.664)* 

-6.250 

(-6.250)* 

-1.904 

(-1.904)* 

20.5 

-19.732 

(-19.732)* 

-17.446 

(-17,446)* 

-13.693 

(-13.693)* 

-8.606 

(-8.606)* 

-2.594 

(-2.594)* 

30.5 

-29.697 

(-29.697)* 

-27.296 

(-27.296)* 

-23.330 

(-23.330)* 

-17.863 

(-17.863)* 

-11.043 

(-11.043)* 

"3‘315* 

(-3.315) 

Note:  Numbers  without  the  parentheses  correspond  to  an  exact  solution 

by  a  numerical  search  technique. 


Table  3.8  Zeros  of  the  Reflection  Coefficient: 


The  Finite  Square  Well 


V  (x)  H 


0 

-E 


o 

0 


;  x  <  -tt/2 

;  -tt/2  <  x  <  tt/2 

;  x  >  tt/2 


Momentum 
\at  which 
\R  is 

E0  zero 

is£ 

2nd 

3rd 

4th 

2.5 

1.22474 

2.54951 

3.67423 

4.74342 

6.5 

1.58114 

3.08221 

4.30116 

12.5 

1.87083 

3.53553 

4.84768 

20.5 

2.12132 

3.93700 

30.5 

2.34521 

4.30116 

X  VALUE 

0.0  0.4  0.6  1  .2  1  .6  2.0  2.4  2.6  3.2  3.6 
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Fig.  3.24.  Results  for  the  Finite  Square  Well  (E0  =  2.5): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 


X  VALUE 
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Fig.  3.25.  Results  for  the  Finite  Square  Well  (EQ  =  6.5): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.26.  Results  for  the  Finite  Square  Well  (E0  =  12.5): 

(a)  Actual  (A)  and  reflectionless  (E)  potential. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 
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Fig.  3.27.  Results  for  the  Finite  Square  Well  (E0  =  20.5): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.28.  Results  for  the  Finite  Square  Well  (E0  =  30.5). 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 
(r)  The  reflection  coefficient. 
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(3.23) 


o 

Again,  as  in  the  case  of  the  harmonic  oscillator  well  and  the  linear 
well,  a  change  in  changes  the  range  a,  too.  Five  different  depths 
were  considered. 


Eigenvalues  and  the  reflection  coefficient  were  calculated  by 
using  the  numerical  methods  of  Appendix  B.  Eigenvalues  and  the  zeros 
of  the  reflection  coefficient  are  summarized  in  Table  3.9  and  Table 
3.10. 


Results  of  the  inverse  calculations  and  the  reflection  coef¬ 
ficients  are  plotted  in  Figure  3.29  to  Figure  3.33. 

The  behaviour  of  the  reflection  coefficient  is  rather  complex 
in  this  case.  The  potential  is  highly  reflecting.  Finally,  with  an 
increase  in  depth  (Eo=80.0  onwards)  the  peaks  of  the  reflection  co¬ 
efficient  seem  to  get  broader. 

The  reflectionless  approximation  seems  to  get  worse  with  an 
increase  in  Eq.  There  is  no  sign  of  convergence  to  the  actual  poten¬ 
tial,  in  much  the  same  way  as  the  square  well. 

3 .2  Confining  Potentials 

As  discussed  in  the  introduction,  the  symmetric  confining  po¬ 
tentials  in  one-dimension  form  a  class  of  their  own  for  the  application 
of  the  inverse  scattering  technique  developed  so  far  because  of  the 
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Table  3.9  Eigenvalues  of  the  Secant-Square  Well 


V  (x) 


0 


>  a 


8 

20 

'•cos  2x 


<  a 


with  E 

o 


\Eigenvalues 

\n  ascend- 

\ng  or- 

E  N^der 

o 

2 

E1  =  -K1 

E2  =  ~K2 

E3  = 

E4  =  A 

E5  =  -k25 

24.0 

(-9.328)* 

(-9.327)+ 

48.0 

(-32.402)* 

(-32.401)+ 

(-14.883)* 

(-14.883)+ 

80.0 

(-64.174)* 

(-64.174)+ 

(-45.160)* 

(-45.160)+ 

(-20.679)* 

(-20.679)+ 

120.0 

(-104.092)* 

(-104.091)+ 

(-84.589)* 

(-84.589)+ 

(-58.196)* 

(-58.196)+ 

(-26.657)* 

(-26.658)+ 

168.0 

(-152.054)* 

(-152.054)+ 

(-132.341)* 

(-132.341)+ 

(-105.228)* 

(-105.228)+ 

(-71.449)* 

(-71.450)+ 

(-32.779)* 

(-32.781)+ 

Table  3.10  Zeros  of  the  Reflection  Coefficient: 


The  Secant-Square  Well 


V  (x) 


0 


''cos  2x 


x  I  _>  a 

with  E 

I  ^  o 

x  <  a 


n.  Momentum  at 
\which  R 
\is  zero 

E  \ 

o  X 

ist 

2nd 

3rd 

-t' 

rt 

I34 

24.0 

1.335 

6.712 

10.486 

14.02 

48.0 

1.273 

6.814 

10.302 

13.443 

80.0 

1.263 

7.23 

10.75 

13.838 

120.0 

1.269 

7.696 

11.328 

14.453 

168.0 

1.2831 

8.162 

11.933 
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Fig.  3.29.  Results  for  the  Secant-Square  Well  (E0  =  24.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.30.  Results  for  the  Secant-Square  Well  (E0  =  48.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.31.  Results  for  the  Secant-Square  Well  (E0  =  80.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The 'reflection  coefficient. 
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Fig.  3.32.  Results  for  the  Secant-Square  Well  (E0  =  120.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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Fig.  3.33.  Results  for  the  Secant-Square  Well  (E0  =  168.0): 

(a)  Actual  (A)  and  reflectionless  (E)  potentials. 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  reflection  coefficient. 
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’ E^-ambiguity ’ .  This  parameter  entered  during  the  translation  of  the 
pure  bound-state  problem  of  a  confining  potential  to  a  suitable  scat¬ 
tering  problem.  To  reiterate  our  plans,  we  will  approximate  a  confining 
potential  V(x),  locally,  by  a  N-bound-state  approximation  V^(x)  ,  where 
[V^(x)  -E^(N)]  is  a  symmetric  reflectionless  potential  supporting  N 

bound-states  at 
2 

— K  =  -E  +  E  ;  n  =  1,2,...,N  .  (3.24) 

n  on  »  >  » 

The  parameter  Eq(N)  will  be  consistently  chosen  to  be 

E  (N)  =  i  (E  +  E  )  (3.25) 

o  2  N  N+l 

where  E^  stands  for  the  n^  bound-state  energy  of  the  confining  poten¬ 
tial. 


3.2.1.  The  Harmonic  Oscillator. 

'  This  is  one  of  the  most  favourite  symmetric  confining  potentials 
in  one  dimension  [45,65].  Choosing  the  parameterization  as. 


V  (x)  =  x‘ 


(3.26) 


the  bound-state  energies  of  this  potential  are  known  to  be 


En+1  "  2n  +  1 


;  n  =  0,1,2,... 


(3.27) 


and  the  bound-state  wavefunctions  are  given  by 


-i/4  2 

00  Hn(x)  e 

n+1  J 2n  n!  n 


(3.28) 


where  H  (x)  is  the  Hermit e  polynomial  given  by 
n 

2  ,n  2 
/1Nnx  ,d.-x 
H  (x)  =  (-1)  e  ( — — )e 

n  dx11 


(3.29) 
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Fig.  3  .34.  Results  for  the  Harmonic  Oscillator: 

(a)  Actual  potential  (A)  and  the  N  = 1  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  corresponding  ’exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.35.  Results  for  the  Harmonic  Oscillator: 

(a)  Actual  potential  (A)  and  the  N  =  2  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  ’exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.36.  Results  for  the  Harmonic  Oscillator: 

(a)  Actual  potential  (A)  and  the  N  = 3  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf .  Section  3.3.2) 


60 


114 


Fig.  3.37.  Results  for  the  Harmonic  Oscillator: 

(a)  Actual  potential  (A)  and  the  N  =  4  reflectionless  approximation  (E)  . 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2)  , 
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Fig.  3.38.  Results  for  the  Harmonic  Oscillator: 

(a)  Actual  potential  (A)  and  the  N  =  5  reflectionless  approximation  (E)  . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 


The  successive  approximations  with  N= 1,2, 3,4  and  5  are  plotted  in 
Figure  3.34(a,b)  to  Figure  3.38(a,b). 


Even  a  casual  look  at  the  five  cases  reveals  an  impressive 
tendency  of  convergence  to  the  actual  potential  with  increase  in  N. 


3.2.2.  The  Linear  Potential. 

We  consider  the  symmetric  linear  potential  given  by  [20,45] 


V (x)  =  lx 


(3.30) 


The  bound-state  energies  are  known  to  be  given  by  the  zeros  of  the 
Airy  functions  [62] 


Ai 1 (-E  )  =  0 
n 


n  = 


1,3,5, 


Ai(-E  )  =  0 

n 


n  = 


2,4,6, 


(3.31) 


The  bound-state  wavef unctions  are  given  by 


ip  (x) 
n 


N. 


Ai(x  -  E  )  ;  x  >  0 

n 


n 


(-1) 


n-1 


N. 


Ai(-x  -  En)  ;  x  <  0 


n 


(3.32) 


The  normalization  constant  is  given  by 


N  =  2 
n 


dx[Ai(x  -  En) ]  =  < 


2E  [ Ai(-E  )]  ;  n  = 1,3,5, .. . 

n  n 


2 [ Ai ’ ( -E  ) ]  ;  n  =  2,4,6, .. . 

n 


(3.33) 


We  display  the  successive  approximations  up  to  N  -  5  in  Figur 
3.39 (a ,b)  to  Figure  3.43(a,b).  There  is  a  strong  evidence  that  the 
local  reflectionless  approximation  converges  to  the  actual  one  as  N 
gets  larger. 


' 


117 


u 


CO 


Fig.  3.39.  Results  for  the  Linear  Potential: 

(a)  Actual  potential  (A)  and  the  N  = 1  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.40.  Results  for  the  Linear  Potential: 

(a)  Actual  potential  (A)  and  the  N  =  2  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.41.  Results  for  the  Linear  Potential: 

(a)  Actual  potential  (A)  and  the  N  =  3  reflectionless  approximation  (E)  . 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf .  Section  3.3.2) 
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Fig.  3.42.  Results  for  the  Linear  Potential: 

(a)  Actual  potentials  (A)  and  the  N  =  4  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  ’exact  reconstruction'.  (cf.  Section  3.3.2) 
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3.2.3.  The  Infinite  Square  Well. 

This  is  again  one  of  the  most  common  examples  in  elementary 
quantum  mechanics.  This  depicts  a  particle  trapped  between  two  im¬ 
penetrable  walls.  In  the  wave  picture,  this  illustrates  the  standing 
waves  [45,61] 


1°  : 

X 

V  (x)  =  - 

OO  • 

y 

X  I 

<  tt/2 
>  tt/2 


(3.34) 


The  bound-state  energies  of  this  potential  are  known  to  be 


E 

n 


2 

n 


n  =  1,2,3 ,.. . 


(3.35) 


The  bound-state  wavefunctions  are 


Vx>  = 


f  r*  -1/2 

TT  cos  nx 


-1/2  . 
tt  sin  nx 


;  n  =  1,3, 

;  n  =  2,4, 


0 


;  |  x  I  <  tt/2 


;  |  x  |  >  tt/2 

(3.36) 


Starting  with  N  =  1,  the  successive  approximations  to  the  infinite 
square  well  up  to  N  =  5  have  been  plotted  in  Figure  3.44(a,b)  to  Figure 
3.48(a,b) . 


Apart  fron  an  oscillatory  stability  at  x  =  0,  the  reflectionless 
approximation  does  not  seem  to  converge  to  the  actual  potential. 

3.2.4.  The  Confining  Poschl-Teller  Potential. 

Let  us  now  consider  the  potential 


k(k  -  1) 

•  2(  17  \ 
sin  (x  —  ) 


X(X-l) 


2  TT. 
COS  (x  -7-) 


V  (x)  = 


(3.37) 


♦ 


- 
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Fig.  3.44.  Results  for  the  Infinite  Square  Well: 

(a)  Actual  potential  (A)  and  the  N  =  1  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.46.  Results  for  the  Infinite  Square  Well: 

(a)  Actual  potential  (A)  and  the  N  =  3  reflectionless  approximation  (E)  . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.47.  Results  for  the  Infinite  Square  Well: 

(a)  Actual  potential  (A)  and  the  N  =  4  reflectionless  approximation  (E)  . 

(b)  Bound-state  wavefunct ions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fie.  3.48.  Results  for  the  Infinite  Square  Well: 

(a)  Actual  potential  (A)  and  the  N  =  5  reflectionless  approximation  (E) 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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For  k  =  A,  this  is  a  symmetric  potential.  This  is,  obviously,  a  perio¬ 
dic  potential.  We,  however,  consider  it  only  between  -  -7-<x<-r.  The 

4  —  —  4 

barriers  put  by  the  singularities  are  impenetrable.  So,  this  suffices 
for  the  confinement  of  a  particle  in  this  potential. 

The  eigenvalues  for  this  potential  are  known  to  be  [61] 

En+1  =  («+A+2n)2  ;  n  =  0,1,2, _  (3.38) 

and  the  unnormalized  wavefunctions  are  given  by 

K  A  i  2 

ipn+1(z)  =  sin  z  cos  z  ^^(-n-jK+A-hijK:  +~2  > sin  z)  (3.39) 

where  z=x-tt/4j  and  is  the  normalization  constant.  We  consider 
this  potential  for  one  specific  value  of  K  and  A,  viz.  k  =  A  =  2.  In 
this  case  eqn.  (3.37)  takes  the  simple  form 

V(x)  =  - ^ —  (3.40) 

cos  2x 

and  eqn.  (3.38)  becomes 

En+1  =  (4  +  2n) 2  ;  n  =  0,1,2,...  .  (3.41) 

Successive  reflectionless  approximations  to  this  confining 
potential  up  to  N  =  5  are  plotted  in  Figure  3.49(a,b)  to  Figure  3.53(a,b). 
Except  for  an  agreement  at  x  =  0,  the  reflectionless  approximation  does 
not  seem  to  exhibit  a  converging  trend.  The  behaviour  is  similar  to 
the  infinite  square  well  which  is,  perhaps,  not  unexpected.  Both  of 
them  have  the  sharp  impenetrable  walls  at  certain  values  of  x . 

3.3.  On  the  Convergence  of  the  Reflectionless  Approximation  to  the 

Actual  Potential. 


Our  investigations  have  left  us  with  a  veriety  of  results.  We 
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Fig.  3.49.  Results  for  the  Confining  Poschl-Teller  Potential: 

(a)  Actual  potential  (A)  and  the  N  = 1  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunction  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.51.  Results  for  the  Confining  Poschl-Teller  Potential: 

(a)  Actual  potential  (A)  and  the  N  =  3  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction'.  (cf.  Section  3.3.2) 
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Fig.  3.53.  Results  for  the  Confining  Poschl-Teller  Potential: 

(a)  Actual  potential  (A)  and  the  N  =  5  reflectionless  approximation  (E) . 

(b)  Bound-state  wavefunctions  of  the  reflectionless  potential. 

(c)  The  corresponding  'exact  reconstruction*.  (cf.  Section  3.3.2) 
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have  observed  cases  of  excellent  agreement  between  the  reflectionless 
approximation  and  the  actual  symmetric  potential.  At  the  same  time, 
however ,  we  have  encountered  instances  of  catastrophic  divergences  from 
the  actual  potential.  Do  we  have,  then,  some  suggestions  regarding 
the  convergence?  We  shall  try  to  answer  this  question,  within  the 
scope  of  the  results  obtained,  in  this  section. 

3.3.1.  Non-conf ining  Potential. 

Things  are  relatively  simple  for  potentials  belonging  to  this 
class  because  of  the  absence  of  the  EQ-ambiguity .  We  have  approximated 
the  actual  potential  by  a  symmetric,  reflectionless  potential  having 

bound-states  at  the  same  energies  as  the  actual  potential.  Exact 

% 

bound-state  energies  were  used  in  the  inverse  reconstruction.  This 
cannot,  then,  be  a  source  of  disagreement  between  the  two.  The  rele¬ 
vant  quantity  to  observe  is  the  reflection  coefficient.  This  was  the 
motivation  for  calculating  the  reflection  coefficient  as  a  function  of 
momentum  or  energy.  The  reflectionless  approximation  is  unique  -  so, 
unless  the  actual  potential  is  reflectionless,  the  two  cannot  be  the 
same.  The  question,  then,  is  that  how  close  the  agreement  is? 

We  saw  one  case  of  fantastic  agreement  -  the  Poschl-Teller 
potential  with  X  =  3.0001.  The  reflection  coefficient  was  equal  to 
unity  at  k = 0  and,  practically,  zero  everywhere  else.  This  suggests 
that  if  the  actual  potential  has  a  reflection  coefficient  which  decays 
extremely  fast  with  increase  in  k,  the  reflectionless  approximation 
should  be  excellent.  This  was  further  vindicated  in  the  case  of  other 
Poschl-Teller  potentials  and  some  of  the  Gaussians,  too. 
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All  other  potentials  except  these  two  had  sharp  'corners'  and 
they  exhibited  pronounced  resonances  in  their  spectrum.  We  have  been 
able  to  obtain  two,  rather  primitive,  suggestions.  First,  if  the  re¬ 
flection  coefficient,  when  plotted  as  a  function  of  momentum, decays 
quickly  then,  despite  the  presence  of  'sharp'  resonances,  the  approxi¬ 
mation  is  tolerable.  Second,  if  the  reflection  coefficient  decays 
very  slowly  and  has  broad  resonances,  the  approximation  is  disappoin¬ 
ting. 

Based  on  these  experiences,  we  can  'qualitatively'  define  two 
classes  of  reflection  coefficients.  If  the  reflection  coefficient 
decays  quickly  and  its  peaks,  if  any,  are  narrow  then  we  will  call 
such  a  reflection  coefficient  to  be  exhibiting  a  favourable  behaviour. 
On  the  other  hand,  the  reflection  coefficient  which  decays  slowly  and 
whose  peaks,  if  any,  are  broad  will  be  said  to  exhibit  an  unfavourable 
behaviour . 

3.3.2.  Confining  Potentials. 

A  convergence  of  the  reflectionless  potential  to  the  actual 
confining  one  at  x  = 0  [47]  does  not  necessarily  guarantee  any  tendency 
towards  convergence  even  in  the  'local'  sense  put  forward  by  Schonfeld 
et.  al.  and  QRT  [45,47].  A  glaring  example  of  such  a  failure  is  the 
Confining  Poschl-Teller  Potential  where  the  approximation  diverges 
despite  its  agreement  at  x  = 0  with  the  actual  potential.  What  is  im¬ 
perative  is  to  obtain  evidences  of  "overall"  convergence  in  the  desired 
region  -  if  possible,  some  criteria  in  support  of  such  a  behaviour. 


Our  investigations  in  the  case  of  non-confining  wells,  which 


lit- 


are  bonafide  members  of  the  class  of  scattering  potentials,  have  promp¬ 
ted  us  to  expect  that  if  the  bound-state  energies  used  in  the  inverse 
method  are  correct  and  if  the  reflection  coefficient  behaves  in  a 
’favourable'  way,  then  the  reflectionless  approximation  is  satisfac¬ 
torily  close  to  the  actual  potential.  Our  numerical  studies  do  not, 
in  totality,  define  the  'favourable'  class  of  reflection  coefficients. 
However,  as  discussed  in  the  earlier  section,  we  have  evidence  for  a 
few  in  this  class.  Undoubtedly,  it  is  beyond  controversy  that  the  two 
quantities  of  interest  in  investigating  the  problem  of  convergence  are 
the  bound-state  energies  and  the  reflection  coefficients. 

•> 

When  we  get  down  to  the  confining  potentials,  there  is  an  in¬ 
herent  artificiality  involved  in  translating  this  purely  bound-state 

problem  into  a  scattering  problem.  As  per  QRT's  prescriptions,  we 

th  th 

introduce  a  zero  of  energy  E^,  half-way  between  the  N  and  (N  +  l) 
bound-states.  We  then  approximate  the  confining  potential  V(x) ,  'lo¬ 
cally',  by  a  potential  V^(x) ,  where  [V^(x)  -Eq]  is  the  reflectionless 
potential  supporting  bound-states  at  the  scaled  'confining  energies' 

(i.e.  bound-state  energies  corresponding  to  the  confining  potential, 

2 

viz.,  E  )  -k  = -E  +  E  .  We  will  call  this  -  the  QRT  approximation. 

*  n  n  o  n 

As  far  as  the  reconstruction  of  confining  potentials  is  concerned,  this 
prescription  is  sufficient  because  the  only  informations  we  have  are 
the  bound-state  energiesof  the  confining  potential.  However,  when  we 
want  to  talk  about  convergence,  we  need  to  carefully  define  the  scat¬ 
tering  counterpart  of  the  confining  potential,  corresponding  to  a  spe¬ 
cific  E  ,  which  we  hope  to  simulate.  That  will,  then,  give  us  some 
o 

information  regarding  the  various  factors  affecting  the  convergence. 


' 


To  put  QRT's  prescription  the  other  way  round  we  scale  down 

the  actual  potential  by  and  try  to  approximate  the  negative  part  of 

the  scaled  potential  by  a  reflectionless  potential  having  bound-states 
2 

at  — E^+E^,  as  shown  in  Figure  3.54(a)  and  Figure  3.54(b).  If  we 
recall  the  properties  of  the  reflectionless  potentials  that  they  are 
always  negative  and  go- to  zero  exponentially  as  x  = ±°°,  we  realize  that 
what  we  can  hope  to  reconstruct,  in  the  best  of  events,  is  the  fol¬ 
lowing  scattering  potential 


V 


scatt . 


V(x)  -  E 

;  |  x  |  <  a 

o 

0 

;  |  x  |  >  a 

(3.42) 

with  the  obvious  equality 


V(±a)  -  E  =  0 
o 


(3.43) 


This  is  the  truncated  potential  shown  in  Figure  3.54(c).  This  is  the 
actual  scattering  counterpart  of  the  confining  potential  V(x)  corres¬ 
ponding  to  a  specific  Eq.  The  successive  approximations  to  the  con¬ 
fining  potential  with  increasing  E^  or  N  then  become  approximations  to 
the  wells  of  the  kind  (3.42)  with  increasing  depth. 


At  the  very  outset,  it  is  clear  that  bound-state  energies  of 
the  scattering  counterpart  are  not  the  same  as  the  scaled  confining 
energies.  The  boundary  conditions  are  to  be  imposed  at  x  = ±a  rather 
than  at  x  =  ±°°.  We  can  expect  the  lower  bound-states  to  be  close  to 
the  assumed  value,  if  the  potential  is  deep  enough,  but  the  upper  ones 
will  be  very  different.  Moreover,  the  discrepancy  between  the  two  is 
highly  dependent  on  the  potential  and  forbids  one  to  draw  any  general 


conclusion. 
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Fig. 


3.54.  Scattering  counterpart  of  a  confining  potential  corresponding 
to  a  specific  EQ. 
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In  the  QRT  reconstruction,  then,  one  tries  to  simulate  a  poten¬ 
tial  with  the  reflectionless  potential  supporting  bound-states  at 
slightly  (hopefully!)  different  energies.  What  is  the  contribution  of 
this  towards  the  convergence  features  of  the  approximation? 

The  obvious  way  to  check  this  is  to  reconstruct  wells  of  the 

kind  (3.42)  with  desired  V(x)  and  E  using  the  true  bound-state  ener- 

o 

gies.  We  call  this  the  exact  reconstruction.  Now  that  we  have  a 
carefully  defined  scattering  counterpart  of  the  confining  potential 
corresponding  to  a  specific  ,  we  can  calculate  its  reflection  coef¬ 
ficient  and  watch  its  behaviour  as  E  increases. 

o 

In  fact,  we  have  already  carried  out  such  a  calculation.  The 
non-confining  wells  in  Section  3.1.3  -  Section  3.1.6  are  precisely  wells 
of  the  type  (3.42)  corresponding  to  the  confining  potentials  in  Section 
3.2.1  -  Section  3.2.4  respectively,  precisely  in  the  same  order.  We 
have  scaled  up  the  results  in  Section  3.1.3  -  Section  3.1.6  by  the  re¬ 
levant  Eq  and  then  displayed  those  ’exact  reconstructions’  with  the 
corresponding  'QRT  reconstructions'  in  Section  3.2.1  -  Section  3.2.4. 

The  reflection  coefficients  for  such  wells  have  already  been  discussed 
in  the  previous  section.  We  will  make  some  remarks  on  them. 

The  Harmonic  Oscillator. 

Looking  at  Figure  3.34(a,c)  to  Figure  3.38(a,c),  we  realize 
that  the  encouraging  convergence  of  the  QRT  reconstruction  has  dis¬ 
appeared  in  the  exact  case.  The  agreement  in  the  exact  case  is,  however, 
tolerable  and  does  not  seem  to  diverge. 


. 
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The  reflection  coefficients  in  this  case  have  already  been 

seen  to  get  favourable  with  an  increase  in  E  . 

o 

The  Linear  Potential. 

A  comparison  between  the  exact  and  the  QRT  reconstruction  can 
be  drawn  from  Figure  3.39(a,c)  to  Figure  3.43(a,c).  The  overwhelming 
convergence  of  the  QRT  reconstruction  has  again  gotten  subdued  in  the 
exact  one.  However,  the  exact  one  shows  some  explicit  tendency  towards 
convergence  unlike  in  the  case  of  the  harmonic  oscillator. 

As  we  have  already  seen  in  Section  3.1.4,  the  reflection  co¬ 
efficient  in  this  case  behaves  more  favourably  than  in  the  case  of  the 
harmonic  oscillator. 

The  Infinite  Square  Well. 

It  is  interesting  to  note  that  the  actual  scattering  counter¬ 
parts  in  this  case  have  not  only  different  eigenvalues  but  they  have 
an  additional  bound-state  as  compared  to  the  numbers  assumed  in  the  QRT 
construction. 

Figure  3.44(a,c)  -  Figure  3.48(a,c)  immediately  suggest  an 
enhanced  diverging  trend  in  the  exact  reconstructions.  The  reflection 
coefficient  in  this  case,  with  its  broad  peaks  and  slow  decay,  has 
already  been  classified  as  unfavourable. 

The  Confining  Posohl-Teller  Potential. 

Figure  3.49(a,c)  -  Figure  3.53(a,c)  confirm  that  the  behaviour 
in  this  case  is  similar  to  the  behaviour  in  the  case  of  the  infinite 
square  well.  With  its  unfavourable  reflection  coefficient,  it  shows  a 


.* 
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definite  divergence  from  the  actual  potential.  The  divergence  is  much 
more  pronounced  in  the  exact  reconstructions. 

This  numerical  exercise  suggests  that  unless  the  reflection 
coefficient  of  the  actual  scattering  counterpart  of  the  confining  po¬ 
tential  with  a  specific  behaves  'favourably T  with  an  increase  in  E^, 
the  convergence  of  the  reflectionless  approximation  to  the  confining 
potential,  even  in  the  local  sense,  is  doubtbul.  The  encouraging  con¬ 
vergence  observed  in  the  QRT  reconstruction  seems  to  be  a  result  of  the 
fact  that  one  uses  the  scaled  confining  energies  there  instead  of  the 
exact  bound-state  energies  of  the  scattering  counterpart  which  one  is 

trying  to  simulate.  This  discrepancy  leads  to  a  scaling  of  the  upper 
• 

% 

levels  in  the  N  bound-state  approximation  relative  to  the  lower  ones. 
The  reconstruction  seems  to  be  quite  sensitive  to  this  scaling. 


4.  DISCUSSION 


As  proposed  in  the  introduction,  the  purpose  of  this  thesis 
was  to  investigate  the  reliability  of  the  reflectionless  approximation 
to  a  symmetric  potential  in  one  dimension  by  applying  it  to  various 
known  potentials  of  this  kind.  Our  results  have  confirmed  that  the 
reliability  of  this  approximation  is  highly  potential-dependent.  This 
is  not  unexpected  as  one  is  trying  to  simulate  a  reflecting  potential 
by  a  reflectionless  one.  Our  calculations  have  unfolded  some  syste- 
matics  in  the  behaviour  of  the  reflection  coefficient  of  the  actual 
potential,  as  a  function  of  momentum  (or  energy),  which  indicate  an 
agreement  (or  a  disagreement)  between  the  two. 

The  question  of  convergence  of  the  reflectionless  approximation 
to  the  actual  potential  is,  again,  a  potential-dependent  query.  An 
agreement  between  the  two  at  a  specific  point  may  be  highly  misleading  - 
it  does  not  necessarily  guarantee  even  a  converging  trend.  What  is  de¬ 
sired  is  an  ’overall’  agreement  between  the  two.  The  only,  obvious, 
way  to  obtain  this  is  to  solve  the  Gel 'fand-Levitan  equation  using  the 
reflection  and  the  transmission  coefficients  of  the  actual  potential  and 
to  show  that,  indeed,  the  difference  between  the  resulting  potential 
and  the  reflectionless  approximation  is  small.  This  is  a  highly  non¬ 
trivial  job  as  one  needs  the  detailed  analytic  structure  of  the  re¬ 
flection  and  the  transmission  coefficients  in  the  complex  k-plane. 

The  problem  of  convergence  ofthe  local,  N-bound  state  approxi¬ 
mation  to  the  confining  potential,  with  increasing  N,  is  still  more 
formidable.  This  approximation  to  the  confining  potential  is  actually 
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an  approximation  to  a  'scattering  counterpart*  with  a  certain  depth. 

This  was  brought  into  light  in  Chapter  3 .  The  scaled  conf ining-energies 
used  in  the  QRT  reconstruction  are  not  the  eigenvalues  of  the  scatter¬ 
ing  counterpart  which  one  is  trying  to  simulate.  In  fact,  they  differ 
considerably.  In  some  cases  even  the  number  of  eigenvalues  of  the 
'scattering  counterpart'  is  different  from  the  number  of  scaled  con- 
fining-energies  -  in  contrast  to  the  assumption  implicit  in  QRT's 
prescription.  We  observed  that  the  effect  of  these  discrepancies  is 
quite  substantial.  So,  there  are  two  crucial  factors  which  govern  the 
convergence  in  this  case:  (i)  a  discrepancy  in  the  number  and  values 
of  the  scaled  conf ining-energies  which  are  used  in the  QRT  reconstruc¬ 
tion;  and  (ii)  the  scattering  counterpart  is,  in  general,  a  reflecting 
potential.  The  eigenvalue  dependence  of  the  reflectionless  potential 
is  rather  involved  and  forbids  us  from  making  any  comment.  Our  inves¬ 
tigations,  however,  show  that  the  successive  'local'  reflectionless 
approximations  to  the  confining  potential,  as  per  QRT's  prescription, 
do  not  show  an  'overall'  convergence  unless  the  reflection  coefficient 
of  the  scattering  counterpart  behaves  in  a  'favourable'  way  with  an 
increase  in  depth. 

From  a  practical  point  of  view,  the  problem  of  convergence  of 
the  local  N  bound-state  reflectionless  approximation  to  a  confining 
potential  with  an  increase  in  N,  is  not  very  demanding.  The  method 
has  computational  limitations,  anyway,  for  large  N.  The  numbers  invol¬ 
ved  become  too  large  or  too  small  to  command  a  faith  in  themselves.  The 
exact  limiting  value  of  N  is  again  dependent  on  the  potenatial.  How— 
ever ,  it  is  enough  to  realize  that  there  is  a  practical  cut-off  to  N. 
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The  important  practical  concerns  are,  then:  (i)  Does  it  show  a  con¬ 
verging  trend  or  a  diverging  trend?;  (ii)  How  fast  is  the  convergence 
or  divergence?.  If  the  reflectionless  approximation  converges  very 
fast  or  diverges  very  slowly  then  the  usefulness  of  this  technique  is 
retained.  On  the  other  hand,  if  it  converges  very  slowly  or  diverges 
very  fast  then  it  is  of  no  real  use. 

_  * 

To  make  a  comment  on  the  overall  performance  of  this  technique, 

one  needs  to  recall  the  kind  of  job  it  was  supposed  to  do.  We  knew 
very  few  energy  levels  of  a  symmetric  potential  in  one  dimension.  That 
was  all.  From  that  bit  of  information,  this  method  was  supposed  to 
generate  the  scattering  potential  for  us.  Even  a  mild  cynicism  towards 
this  technique  will  accede  to  the  fact  that  the  input  information  is 
too  meagre  to  expect  a  grand  result.  Within  this  limitation,  the  per¬ 
formance  of  this  technique  has  been  satisfactory.  For  a  small  number 
of  eigenvalues,  which  is  its  actual  domain  of  application,  it  simulates 
the  rough  'qualitative'  features  of  the  potential  quite  well.  Even  in 
the  worst  cases  which  we  encountered,  it  showed  at  least  this  much 
promise.  Any  ’ quantitative' expectation  from  this  technique  is,  perhaps, 
highly  ambitious.  This  is,  indeed,  beyond  the  scope  of  any  such  naive 
technique . 


. 
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APPENDIX  A.  SOME  IMPORTANT  LEMMAS 


Lemma  1.  The  determinant  of  the  matrix 


Kn+Kv 


is  given  by 


=  det 


K  +K 

n  v 


N  2  N  2  N  2 

M  n  (k ,-k.)  n  (k0-k.) —  n  (k.t  --0 

N  1  .1  ,1>1  1  J  r>2  2  l  1>N-1  N~X  3 


n  OtH) 


N 


N 


^  J  N  r\  2  *1  2 

n  (k.+k.)  n,  (k,+k.)  ....  n  (<„.+<.) 

j>X  1  J  j>2  2  J  j  >N-1  N_1  J 


with  n = V  = 1,2 , . . . . ,N. 


Proof.  Let  us  define 


A  =  det 

p  "  K  +K 

*  n  v 


;  n  =  v  =  1, . .  .  ,p 


(A.  1) 


We  will  then  prove  this  lemma  by  induction  on  p.  For  p  =  l. 


A-  =  — t - 

1  1  k^+k^ 


2k. 


(A. 2) 


For  p  =  2 , 


A2  = 


<2«i 


k1HC2 


K2+K2 


(A. 3) 


Subtracting  the  first  row  from  the  other 

1 


A„  = 


(KrK2) 


2  CK1+K1)  (|C1'h<2') 


k2*1 


K2+K2 


(A. 4) 


Subtracting  the  first  column  from  the  others 

2  2 

i  i  «r<2>2  r 2  i  i 

A2  =  (2r)(2T->  11  4 - 

Z  ZK1  ZK2  (k,+k0)  ln=l  n  J  \  2 


1  2 


n  (k,+k  ) 
j>i  1  j 


(A. 5) 
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If  we  do  exactly  the  same  thing  for  p  =  3 ,  we  get 


3  2  3  2 

n  (k,-k.)  n  (k  -<.) 

j>l  2  j  >2  1  J 

3  2  3  2 

n  (k,+k.k  n  (k9+k.)z 

j>l  J  j>2  Z  J 


(A. 6) 


We  assume  the  result  for  p=N-l,  i.e. 


N-l 

n  ( 

n=l 


N-l  2n-x  2 

n  (<,-<.)  n  (k9-«:.)  . 

j>l  J  j>2  z  J 

N-l  •  N-l 

n  (k  +k.)  n  (k  +k.)z 
j>l  3  j>2  z  3 


N-l 

.  n 

J>N-2 

N-l 

n 


j  >N-2 


(kn-2-V 


(A. 7) 


Finally,  we  prove  it  for  p=N 


V*! 


V*1  Kl+K2 


K.+K  k2+K2 


V*2 


<1^N 


<2+<N 


VSi 


Subtracting  the  first  row  from  all  the  others,  we  get 


N 

n  (k  -k.) 

i>i  1  j 

N 

2k  n  (k  -h<.) 

ij>i  j 


i 

i 


i 

i 


i 

i 


K2-H<1 


k2+k:2 


i<2"H<N 


V*1 


KN+K2 


V*H 


(A. 8) 


Subtracting  the  first  column  from  all  the  others,  we  have 


I  • 
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Lemma 


Proof . 


N  2 

n  (<,-<. r 

j>i  1  J 

N  f  2 

2Ki.n(Ki+Kj) 

j>1  J 


0 

1 


0 

1 


k2+<2 


K2^N 


KN+K2 


kn'hcn 


(A. 9) 


N  2 

n  (k  -k.)z 

i>i  1  j 

N 

2k  n  (k  +k.) 

j>i  1  3 


Vi(n,v 


=  2,....,N) 


r  n 

;  ",  (2^> 

11=1  n 


N  N  N 

n  (k -o  n  (k  -k.)  ....  n  (km  -k.) 

3>1  1  3  j>2  2  3 _ j>N-l  "-1  3  , 

N  N  2  N  2  ^  ; 

n  (k  -he.)  n  (k_+k . ) z _  n  ( k  -h< . ) ^ 

j>l  1  3  j>2  2  3  j>N-l  N-X  3 


QED  . 


1.  It  can  be  proved  that 


N 

n 

n=m=l 

n^ 


K  -K 
n  m 


K  +K 
n  m 


N  2  N  2 

II  (K  -K.)  II  (K  -K.)  .. 
j>l  2  j  >2  J 


N 

n  (k„  -k.) 


j  >N-1 


N-l  j 


N  2  N  ? 

n  (Kt+k.)  n  (k„+k.)  .. 
j>l  1  3  j>2  2  3 


N 

'  11  (KN-1+Ki) 

j  >N-1  J 


We  define 


n  =  n 

(p)  n,m=l 
n?^m 


K  — < 

n  m 


K  +K 
n  m 


(A. 11) 


We  will  then  prove  this  lemma  by  induction  on  p.  For  p  =  l,  the  equality 
of  the  two  products  in  the  statement  is  obvious.  For  p = 2; 
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2 

n 

n,m=l 

n^m 


K  — (C 
n  m 

k1-<2l 

l^l 

k  -He 
n  m 

1  Kj+Kj  | 

Ik2+K1I 

n  (k  -k.) 


 j>i 


i  j 


n  (k  +k.) 
j>l  1  3 


For  p  =  3; 


3 

n 

n  ,m=l 
n^m 


K  — K 
n  m 


K  +K 
n  m 


K1"K2 


x< 


K1*K3 


K1+K2I  l,el'tK3 


/• 

K3-Kil 

1 <3-K2 1 ' 

K.-HC.  I 

1 <„+K^ 1 

L 

3  l1 

1  3  2 1 

y 

K2-k1 

K2+Kl 


K2  K3 

k2+<3 


n  (k.-k.)  n  (k  -k.) 
=  j>l  1  J  ,1>2  2  J 


n  (k.-hc.)  n  (k  +k.) 
j>l  1  J  j>2  2  J 


2 


2 


We  assume  this  result  for  p=N-l,  i.e.. 
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j  >N-2 


(KN-2+Kj) 


Finally,  we  prove  it  for  p  =  N 
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k  -He 
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k2+KN 
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VKi!  Iv^L.JWi 


Wl 


f  N-l 

11 1 
[n,m=l 

n#n 
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(KiV2 
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(A. 12) 


(A. 13) 


(A. 14) 


a  v  a  2 
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2  N  ^  N  2  ^  2 

(k  -k.)  n  (k  -K.)  —  n  (k  -k.)  n  (k  ,-k.) 

1  1  3  ,j>2  2  J  j>N-2  N~2  3  j>N-l  N 1  J 


N 


N 


N 


(Kt+k.)^  n  (k,+k.)2 _  n  (<  0+< , ) 2  n  (k  ,+k.)2 

j>l  1  3  j>2  2  J  j>N_2  N-2  3  j>N_1  N-l  3 


QED  . 


(A. 15) 


APPENDIX  B.  NUMERICAL  METHODS  FOR  CALCULATING  THE  REFLECTION  COEFFICIENT 
AND  THE  EIGENVALUES  OF  A  SYMMETRIC  POTENTIAL  IN  ONE  DIMENSION. 


Eigenvalues 

Two  methods  were  used  to  numerically  calculate  the  eigenvalues. 
The  first  one  is  fairly  general  and  does  not  need  the  specific  assump¬ 
tion  of  symmetry  of  the  potential.  We  discuss  that  first 


(i)  The  Matrix  Method:  The  Schrodinger  equation  for  a  poten¬ 
tial  V(x)  in  one  dimension  is 


-Vl l)  =  -Elp 
dx 


•  _oo  <  X  <  00 


(B.l) 


where  E  is  the  energy  of  the  particle.  We  have  to  find  those  solutions 
which  satisfy  the  boundary  conditions 

1K+00)  =  =  0  .  (B  .2) 


For  numerical  computations  we  need  to  consider  a  finite  interval 
a  <_  x  _<  b.  We  assume  the  boundary  conditions  to  be  ip(a)  =  ip(b)  =  0. 
Physically,  this  means  that  we  are  putting  infinite  impenetrable  walls 
at  x  =  a  and  x  =  b.  We  divide  this  interval  a  <_  x  _<  b  into  (n  + 1)  fine 
intervals  of  width  fh',  i.e. 


h  = 


b-a 

n+1 


where 


a  +  (n+l)h  =  b 


(B  .3) 


(B.4) 


If  we  replace  the  Schrodinger  equation  (B.l)  by  its  corresponding  dif¬ 
ference  equation,  we  get 
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Using  the  boundary  conditions  at  x  =  a  and  x  =  b,  we  get  the  following 
set  of  equations  from  (B.5) 


(B.6) 

where  we  have  multiplied  throughout  by  a  minus  sign.  ..So,  the  eigen¬ 
value  problem  for  the  Schrodinger  equation  is  now  reduced  to  this 
algebraic  eigenvalue  problem  for  the  tridiagonal  matrix  in  eqn.  (B.6) 
The  lowest  negative  eigenvalues  give  us  the  required  bound-state 
energies  for  the  Schrodinger  equation.  These  were  calculated  using 
the  EQRT1S  subroutine  from  the  IMSL  Fortran  Library. 

In  practice,  the  computation  goes  in  two  steps.  First,  one 
picks  up  a  particular  pair  Ta'  &  'b'  and  repeats  the  calculation  of 
the  eigenvalues  by  increasing  the  number  of  intervals  (n  + 1)  till  one 
gets  a  value  stabilized  to  a  desired  accuracy.  Then  one  increases  'a 
and  'b'  and  repeats  the  same  procedure.  When  the  stabilized  values 
for  two  consecutive  'af  &  'b*  agree  within  a  desired  tolerance,  one 
stops. 

The  advantages  of  this  technique  are  its  simplicity  and  the 
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fact  that  it  can  be  used  for  any  kind  of  potential.  However,  the 
process  of  stabilizing  the  eigenvalue,  even  to  a  few  places,  is  very 
tedious.  It  serves  as  a  method  to  quickly  locate  the  positions  of  all 
the  levels  at  the  same  time. 


(ii)  Solving  the  Eigenvalue  Equation  for  a  Symmetric  Potential 
Numerically:  Let  us  consider  a  symmetric  potential  of  the  form 


0 


v(x)  =  mx) 


o 


;  x  -a  (Region  I) 

;  -a  _<  x  <_  a  (Region  II) 

;  x  >  a  (Region  III) 


(B.7) 


and  V(x)  =V(-x).  Also,  it  is  obvious  that  V(±a)  =0.  We  can,  then, 
immediately  write  down  the  solutions  of  the  Schrodinger  equation, 
ip  and  ij;  in  region  I,  II  and  III,  respectively  (for  E  <  0)  : 


where 


ipT  =  Ae 


kx 


ipT-  =  Bip  +  Cip 
rII  Te  To 


W  De 


-kx 


(B  .8) 


E  =  k 


(B  .9) 


Since  the  potential  is  symmetric  in  Region  II,  the  two  linearly 
independent  solutions  in  this  region  are  the  odd  (^Q)  ant^  even  (t^) 
solutions.  If  we  now  impose  the  continuity  conditions  at  x = -CL  and 
use  the  symmetry  properties  of  3.nd  ,  we  immediately  get  that  the 
even  eigenvalues  are  the  solutions  of  the  equation 


ip '  (x  =  +a)  +  kip  (x  =  +a)  =  0 
e  e 


(B.10) 
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and  the  odd  eigenvalues  are  the  solutions  of  the  equation 

^(x=+a)  +  ki|>  (x«+a)  =  0  (B.ll) 

These  two  equations  were  solved  numerically  to  find  out  the  values  of 

the  bound-state  energies.  Values  of  i b  .  \ b  ,1/’  and  ilj '  at  x  = +ct  were 

e  ro  e  o 

needed  in  this  process.  They  were  calculated  using  the  DVERK  subroutine 
from  the  IMSL  Fortran  Library  with  the  obvious  conditions 

^'(0)  =  0 

Ip  (0)  =  0  .  (B .  12) 

o 

This  subroutine  uses  the  Runge-Kutta  method  to  do  this  job. 


This  technique  is  much  more  accurate  and  convenient.  Note: 

*•  * 

For  potentials  with  infinite  tails  like  the  Gaussian,  the  same  method 
was  used.  Some  value  of  a  was  assumed  to  start  with.  Then  it  was 
gradually  increased  and  the  calculations  were  repeated  each  time.  One 
stopped  when  two  consecutive  calculations  for  the  eigenvalues  agreed 
within  a  desired  accuracy. 


Reflection  Coefficient 

We  concentrate  on  a  well  of  the  type  (B.7).  We  are  now  con- 

2 

cerned  with  the  continuum  solutions  (E  = k  )  >0.  Solutions  in  the  three 

classified  regions  are 

,  .  ikx  .  -ikx 

\p  =  A+e  +  A_e 

^II  =  V'e(x)  +  BA(X) 


*  =  C  eikx 

III  + 


(B .  13 ) 


The  other  solution  e  in  is  absent  because  there  is  no  incoming 

wave  in  region  III.  ip  (x)  and  ^  (x)  are  again  the  even  and  odd  solu- 
tions  for  E  >  0  in  region  I.I.  Imposing  the  continuity  conditions  at  • 
x  =  ±a  and  using  the  symmetry  properties  of  the  solutions  ipQ(x)  and 
ipQ(x)  one  immediately  gets  the  reflection  coefficient 


R  = 


rip '  (x=-Hx)ipt  (x=+a) 
e  o 


-  k 


+  ki^  (x=+a)ip  (x=+a) 
e  o 


w2+ 


ip '  (x=+a)  ip  ’  (x=+a) 
e  o 


(B .14) 


+  ki^  (x=+a)ii  (x=+a) 
e  o 


where  W  is  the  Wronskian 

W  =  ip  (x)ip’  (x)  -  ip  (x)^’(x) 
o  e  e  o 


(B .15) 


The  Wronskian  can  be  easily  evaluated  at  x  =  0.  For  a  numerical  calcu¬ 
lation  of  R  as  a  function  of  momentum  k,  one  needs  the  values  of  ip^, 
ip ip  ,  ip*  at  x=+a.  These  were  again  calculated  using  the  DVERK  sub¬ 
routine  which  uses  the  Runge-Kutta  method. 

Note:  For  potentials  with  infinite  tails,  one  starts  with  an  assumed 

value  for  a.  It  is  then  gradually  increased  and  the  calculations  are 
repeated.  Once  the  value  of  R  stabilizes  to  a  certain  desired  accuracy 


one  stops. 


» 


